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14.  ABSTRACT  The  overall  objective  of  this  project  is  to  test  the  hypothesis  that  HOXC9  expression  levels  have  a  causal  role  in 
determining  the  differentiation  states  of  neuroblastoma  tumors,  with  higher  levels  of  HOXC9  promoting  differentiation.  At  the  cellular 
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grant,  significant  progress  has  been  made  in  the  proposed  studies,  which  is  described  in  detail  in  the  Overall  Project  Summary  section 
below.  Briefly,  we  have  2  publications  reporting  our  new  findings  on  the  molecular  mechanism  by  which  HOXC9  induces  neuronal 
differentiation  of  neuroblastoma  cells.  We  show  that  HOXC9  directly  regulate  a  large  number  of  genes  involved  in  neuronal  differentiation, 
upregulating  neuronal  genes  and  downregulating  cell  cycle  and  DNA  repair  genes.  For  Aim  1,  we  have  completed  the  experiments 
proposed  in  MYCN  mouse  tumor  development  studies.  For  Aim  2,  we  have  developed  a  culture  system  for  enriching  and  long-term 
propagating  neuroblastoma  stem  cell  and  shown  that  high  Hoxc9  expression  is  essential  for  their  differentiation  induced  by  RA.  For  Aim  3, 
we  have  identified  3  modes  of  epigenetic  regulation  of  transcription  by  HOXC9  that  involve  modulation  of  H3K4me3  and  H3K27me3 
levels,  thereby  providing  a  molecular  basis  for  the  role  of  KDM5B  and  KDM6B  in  the  control  of  neuroblastoma  differentiation.  These 
findings  significantly  advance  our  molecular  understanding  of  neuroblastoma  differentiation  and  suggest  new  targets  for  neuroblastoma 
therapy. 
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INTRODUCTION: 

Neuroblastoma  is  a  common  pediatric  cancer  of  the  sympathetic  nervous  system  [1],  It  is  a  heterogeneous  group 
of  tumors,  ranging  from  tumors  composed  predominantly  of  poorly  differentiated  neuroblasts  to  those 
consisting  largely  of  differentiated  neurons.  Patients  with  poorly  differentiated  neuroblastoma  have  a 
significantly  poorer  prognosis  than  those  with  differentiated  neuroblastoma.  The  molecular  mechanisms 
underpinning  neuroblastoma  heterogeneity  are  largely  unknown  [2-5].  Our  hypothesis  is  that  HOXC9 
expression  levels  have  a  causal  role  in  determining  the  differentiation  states  of  neuroblastoma  tumors,  with 
higher  levels  of  HOXC9  promoting  differentiation  [6,  7].  At  the  cellular  level,  HOXC9  promotes  the 
differentiation  and  represses  the  self-renewal  of  neuroblastoma  stem  cells.  At  the  molecular  level,  HOXC9 
activates  the  H3K4  demethylase  KDM5B  and  the  H3K27  demethylase  KDM6B  for  global  control  of  its 
differentiation  program.  The  study  supported  by  this  award  has  three  specific  aims:  1)  to  investigate  the  role  of 
HOXC9  in  neuroblastoma  development;  2)  to  investigate  the  role  of  HOXC9  in  neuroblastoma  stem  cell 
differentiation;  and  3)  to  investigate  the  molecular  mechanism  for  global  control  of  HOXC9-induced 
differentiation.  Identification  and  characterization  of  neuroblastoma  stem  cells  will  advance  our  understanding 
of  neuroblastoma  heterogeneity  and  provide  a  key  cellular  target  for  earlier  detection  of  neuroblastoma,  for 
better  prediction  of  clinical  outcomes,  and  for  cancer  stem  cell-based  drug  discovery.  A  molecular 
understanding  of  HOXC9-induced  differentiation  will  open  new  avenues  for  the  development  of  more  effective 
differentiation-based  neuroblastoma  therapies. 

1.  KEYWORDS: 

Cancer  stem  cells,  Chromatin  immunoprecipitation  and  sequencing  (ChIP-seq),  Epigenetics,  Histone  H3 
methylation,  Histone  H3  demethylases  (KDMs),  HOXC9,  Mouse  model  of  neuroblastoma  (MYCN  mice), 
neuroblastoma  differentiation;  Retinoic  acid  (RA) 

2.  OVERALL  PROJECT  SUMMARY: 

Task  1.  To  investigate  the  role  of  HOXC9  in  neuroblastoma  development  (months  3-26) 

The  goal  of  Task  1  is  to  test  the  hypothesis  that  HOXC9  expression  levels  have  a  causal  role  in  determining  the 
differentiation  states  of  neuroblastoma  tumors.  We  will  examine  the  effects  of  Hoxc9  deficiency  and 
heterozygosity  on  the  differentiation  states  of  neuroblastoma  tumors  developed  in  MYCN  mice,  an  animal 
model  of  the  human  disease.  As  differentiation  states  affect  tumor  development,  we  will  also  examine  the 
effects  of  Hoxc9  deficiency  and  heterozygosity  on  neuroblastoma  initiation  and  progression  in  MYCN  mice.  A 
total  of  492  mice  will  be  used  for  breeding  (n  =  144)  and  for  the  proposed  hyperplasia  (n  =  108)  and  tumor 
development  (n  =  240)  studies. 

The  Task  1  experiments  outlined  in  the  approved  SOW  for  the  second  budget  year  include  tumor 
development  studies  and  microarray  analysis. 

Task  1.2.  Tumor  development  studies  (months  13-24):  40  mice  for  each  group  (6  groups,  n=240)  will  be 

monitored  for  tumor  development.  Animals  that  are  moribund  will  be  euthanized  and  subjected  to 
full  autopsy.  Tumors  will  be  removed,  measured  for  weight  and  sectioned  for  H&E  and 
immunofluorescence  staining  for  nestin,  Phox2B  and  TH.  The  stained  sections  will  be  examined 
with  regular  and  confocal  fluorescent  microscopes  and  positive  cells  will  be  quantified  from 
randomly  selected  fields.  We  will  also  quantify  tumor  cells  expressing  nestin,  Phox2B  and/or  TH  by 
immunoflow  cytometry. 

We  have  finished  the  collection,  section  and  H&E  staining  of  tumors  developed  in  Hoxc9+/+  and  Hoxc9+/'  mice 
with  MYCN.  We  have  also  completed  immunofluorescence  staining  and  quantification  of  tumor  cells  that 
express  nestin  (stem  cell  marker),  Phox2B  (progenitor  cell  marker)  and/or  TH  (differentiation  marker).  We  are 
currently  in  the  process  of  analyzing  the  data  to  determine  the  effects  of  Hoxc9  heterozygosity  on 
neuroblastoma  development  and  differentiation. 

We  have  not  been  able  to  obtain  tumors  from  Hoxc9“  MYCN  mice  because  most  of  them  died  within  one 
month  (~4  weeks)  after  birth.  Necropsy  revealed  no  overt  tumor  growth.  In  general,  MYCN  mice  develop 
neuroblastoma  tumors  at  a  median  age  of  ~14  weeks  [8].  We  are  currently  investigating  the  cause  of  the  early 
death  of  Hoxc9 '  '  MYCN  mice. 
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Task  1.3.  Microarray  analysis  of  tumor  differentiation  states  (months  20-26):  we  will  conduct  gene  profiling  of 
Hoxc9'  ',  Hoxc9+/',  and  Hoxc9+/+  tumors  by  microarray.  Five  tumor  samples  per  genotype  will  be 
analyzed  (total  15  microarrays).  We  will  then  analyze  the  microarray  data  using  GO  terms  analysis, 
IPA  and  GSEA  to  gain  a  molecular  view  of  tumor  differentiation  states,  as  well  as  genes  regulated 
by  Hoxc9. 

We  have  collected  Hoxc9+/+  and  Hoxc9+/'  tumors  samples  and  finished  the  preparation  of  RNA  samples  for 
microarray  analysis.  We  are  now  conducting  microarray  assays.  We  anticipate  the  completion  of  microarray 
data  analysis  in  the  next  two  months. 
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Fig.  1.  Isolation  and  functional  characterization 
of  MYCN  neuroblastoma  sphere-forming  cells. 
A-B,  Serial  sphere  formation  assays  of  neuro¬ 
blastoma  cells  from  Hoxc9+/+  MYCN  mice. 
Primary  spheres  were  individually  collected, 
dissociated  into  single  cells  and  plated  into  96- 
well  plates  with  an  average  of  one  cell  per  well. 
C,  Tumorigenic  assays  of  primary  and  sphere¬ 
forming  neuroblastoma  cells. 


Task  2.  To  investigate 
the  role  of  HOXC9  in 
neuroblastoma  stem  cell 
differentiation  (months  16-36) 

The  goal  of  Task  2  is  to  test  the 
hypothesis  that  HOXC9 
expression  levels  control  the 
balance  between  the 
differentiation  and  self-renewal 
of  neuroblastoma  stem  cells.  We 
will  examine  the  effects  of 
Hoxc9  deficiency  and 
heterozygosity  on  the  abilities  of 
nestin'"  neuroblastoma  stem  cells 
to  self-renew  and  to  differentiate 
in  culture,  and  to  generate 
tumors  that  recapitulate  the 
heterogeneity  of  the  original 
tumor  in  syngeneic  mice.  The 
Task  2  experiments  outlined  in 
the  approved  SOW  for  the 
second  budget  year  include 
tumor  cell  sphere  formation  and 
differentiation  assays. 

Task  2.1.  Self-renewal 
sphere  or  soft  agar  colony 
formation  assays  (months  16- 

21):  Tumor  samples  from  Task 
1  will  be  used  to  isolate  nestin- 
GFP+  neuroblastoma  cells  (5 
tumor  samples  per  genotype,  3 
genotypes,  Hoxc9+/+,  Hoxc9+/', 
Hoxc9 '  '  mice  with  MYCN  and 
nestin-GFP).  The  freshly 
isolated  tumor  cells  will  be 
plated  at  1  cell/well  (one  96- 
well  plate  per  tumor  sample)  for 


sphere-formation  assay  or  at  500  cells/well  (one  6- well  plate  per  tumor  sample)  for  colony  formation  assay. 
This  will  be  followed  by  secondary  and  tertiary  sphere  and/or  colony  formation  assays  with  5  randomly 
selected  primary/secondary  spheres  or  colonies  for  each  tumor  sample  (5  tumor  samples  per  genotype). 
Self-renewal  will  be  quantified  as  the  number  of  spheres  or  colonies  generated  from  each  previous  sphere  or 
colony. 
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We  have  performed  all  the  assays  proposed  under  this  task  for  Hoxc9+/+  and  Hoxc9+/‘  MYCN  tumors  (as  stated 
above,  we  have  not  been  able  to  obtain  Hoxc9 '  '  MYCN  tumors). 

We  have  developed  a  cell  culture  system  that  allowed  us  to  propagate  neuroblastoma  sphere-forming  cells 
(Fig.  1A).  Through  serial  sphere  formation  assays  with  single  cells,  we  were  able  to  demonstrate  that  these 
sphere-forming  cells  possess  self-renewal  potential  (Fig.  IB),  an  essential  property  of  cancer  stem  cells. 
Compared  to  their  parental  primary  tumor  cells,  these  sphere-forming  cells  showed  a  marked  (>  100-fold) 
increase  in  the  capacity  to  induce  tumors  in  syngeneic  wild-type  129/J  mice  (Fig.  1C).  Together,  these  data 
indicate  that  we  have  isolated  a  minor  population  of  mouse  neuroblastoma  cells  with  cancer  stem  cells 
properties. 

Importantly,  neuroblastoma  tumors  from  Hoxc9+/'  mice  appeared  to  contain  a  higher  number  of  sphere- 
forming  cells  (~l-5%  of  primary  tumor  cells)  in  comparison  to  Hoxc9+/+  MYCN  tumors  (-0.01%),  suggesting 
that  a  decrease  in  Hoxc9  expression  enhances  the  production  and/or  maintenance  of  neuroblastoma  stem  cells, 
probably  as  a  result  of  reduced  differentiation.  This  finding  is  consistent  with  our  model  that  HOXC9  has  an 
important  role  in  promoting  neuroblastoma  differentiation  [6,  9]. 

Although  Hoxc9'/‘  MYCN  tumors  are  not  available  at  the  moment,  we  have  developed  an  alternative 
approach  to  investigating  the  effect  of  Hoxc9  deficiency  on  the  self-renewal  and  differentiation  of 
neuroblastoma  stem  cells.  We  have  begun  the  process  of  generating  Hoxc9  knockout  neuroblastoma  sphere¬ 
forming  cells  using  the  RNA-guided  CRISPR-Cas9  nuclease  system  [10].  Once  generated,  the  Hoxc9 '  ' 
neuroblastoma  sphere-forming  cells  will  be  examined  for  their  ability  to  undergo  self-renewal  and 
differentiation,  in  comparison  with  their  parental  Hoxc9+/+  sphere-forming  cells. 
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Fig.  2.  Activation  of  the  serine-glycine  synthesis  pathway  in  sphere-forming  cells.  A.  GO  analysis 
of  genes  upregulated  in  sphere-forming  cells  (enrichment  fold  >  2.0,  p  <  0.01).  B,  Serine-glycine 
synthesis  pathway  with  the  indicated  fold  changes  in  mRNA  expression  of  the  enzyme  genes 
determined  by  microarray.  3-PG,  3-phosphoglycerate;  3-PHP,  3-phosphohydroxy-pyruvate; 
3-PSer,  3-phosphoserine;  OAA,  oxaloacetate;  Ptd-Ser,  phosphatidylserine. 


To  characterize  sphere¬ 
forming  neuroblastoma  cells 
at  the  molecular  level,  we 
performed  microarray 
profiling  of  these  cells  from 
3  Hoxc9+/+  MYCN  tumors  in 
comparison  with  their 
parental  primary  cells.  Gene 
Ontology  (GO)  analysis 
revealed  that  among  the 
genes  upregulated  in  sphere¬ 
forming  cells,  those  within 
the  serine-glycine  synthesis 
pathway  were  significantly 
enriched  (Fig.  2A),  including 
phosphoglycerate 
dehydrogenase  (PHGDH), 
phosphoserine 

aminotransferase  1  (PSAT1), 
phosphoserine  phosphatase 
(PSPH),  and  serine 
hydroxymethyltransferase  2 
(SHMT2)  (Fig.  2B).  It  is 
known  that  increased 


activation  of  the  serine-glycine  biosynthetic  pathway  is  an  important  part  of  cancer  metabolism  [11,  12].  In 
addition  to  generating  serine  and  glycine  for  the  biosynthesis  of  proteins,  purines  (via  one-carbon  units),  and 
lipids  (via  phosphatidylserine),  this  pathway  produces  equimolar  amounts  of  reduced  nicotinamide  adenine 
dinucleotide  (NADH),  a-ketoglutarate  (a-KG),  and  5,10-MTHF  (Fig.  2B).  These  metabolites  have  critical  roles 
in  the  control  of  cellular  metabolism  for  cell  proliferation  and  survival:  NADH  participates  in  ATP  production 
and  redox  regulation  [13];  a-KG  supplies  carbon  to  the  tricarboxylic  acid  (TCA)  cycle  for  the  generation  of 
many  essential  biosynthetic  precursors  [14];  and  5, 10-methylenetetrahydro folate  (5,10-MTHF)  is  a  coenzyme 
for  the  only  cellular  pathway  of  de  novo  thymidylate  biosynthesis  catalyzed  by  thymidylate  synthase  and  a 
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major  source  of  one-carbon  units  for  purine  synthesis  [11,  15-17].  Our  findings  provide  the  first  line  of  direct 
evidence  for  increased  activation  of  this  biosynthetic  pathway  in  cancer  stem  cells. 


Task  2.2.  Differentiation  assays  (months  21-26): 

a.  For  spontaneous  differentiation  assays,  clonally  derived  secondary  and  tertiary  spheres  or  colonies  (n 
=  5  per  tumor  sample,  5  tumor  sample  per  genotype,  3  genotypes,  Hoxc9+/+,  Hoxc9+/-,  Hoxc9-/- 
mice  with  MYCN  and  nestin-GFP)  will  be  cultured  at  1  sphere  or  colony  per  coverslip.  The  spheres 
or  colonies  will  be  cultured  for  8  days  and  then  stained  for  Phox2B  or  TH  (neuronal  markers)  and 
BLBP  or  SI 00  (both  are  glial  markers).  Percentages  of  spheres  or  colonies  containing  both  neurons 
and  glial  cells  or  only  neurons  or  glial  cells  will  be  determined. 

b.  For  induced  differentiation  assays,  clonally  derived  secondary  and  tertiary  spheres  or  colonies  will 
be  cultured  at  1  sphere  or  colony  per  coverslip  in  the  presence  of  RA  (neuronal  differentiation). 

After  7  days,  cells  will  be  stained  for  Phox2B  and  TH.  The  percentage  of  neuron-only 
spheres/colonies  will  be  determined.  In  parallel  experiments,  we  will  treat  the  cells  with  neuregulin 
(glial  differentiation).  After  7  days,  cells  will  be  stained  for  BLBP  and  SI  00.  The  percentage  of  glia- 
only  spheres/colonies  will  be  determined  (months  21-26). 


Sphere_RA 


Hoxc9+/+_DMSO 


Fig.  3.  Hoxc9+/-  spheres  are  resistant  to  RA-induced  differentiation. 

A,  Morphological  examination  of  Hoxc9+/+  and  Hoxc9+/-  spheres  treated 
with  DMSO  or  10  ijM  RA  for  7  days.  B,  Quantitative  RT-PCR  analysis  of 
mRNA  expression  of  the  indicated  genes  in  Hoxc9+/+  and  Hoxc9+/- 
spheres  following  treatment  with  DMSO  or  1 0  /vM  RA  for  7  days.  Data 
were  analyzed  by  two-tailed  Student’s  t-test  with  P  values  indicated. 


We  originally  proposed  to  use  Phox2B  and 
TH  as  neuronal  markers  for  assessing 
neuronal  differentiation  of  neuroblastoma 
cells.  However,  our  recent  studies  have 
shown  that  upregulation  of  the  neuronal 
genes  GFRA3  and  NEFM  correlates  better 
with  neuroblastoma  neuronal  differentiation 
[6,  9]  (see  also  Appendix  1).  GFRA3 
encodes  glial  cell  line-derived  neurotrophic 
factor  (GDNF)  family  receptor  alpha  3 
(GFRa3)  and  has  a  critical  role  in 
embryonic  development  of  the  sympathetic 
nervous  system,  promoting  the 
differentiation  and  axonal  outgrowth  of 
sympathetic  neurons  [18].  NEFM 
(neurofilament,  medium  polypeptide 
150kDa)  is  a  common  marker  for 
differentiated  neurons.  Therefore,  we  have 
been  using  Gfra3  and  Nefm  as  markers  for 
neuronal  differentiation  in  both 
spontaneous  and  induced  differentiation 
studies. 

We  have  completed  the  spontaneous 
differentiation  study.  We  found  no 
significant  levels  of  spontaneous 
differentiation  for  both  Hoxc9+/+  and 
Hoxc9+/'  spheres  following  8-day  culture, 
as  determined  by  immunofluorescent 
staining  for  neuronal  (Gfra3  and  Nefm) 
and  glial  (BLBP  and  SI 00)  markers,  and  by 
quantitative  RT-PCR  analysis  of  mRNA 
expression  of  these  markers  (data  not 
shown). 


We  have  also  completed  RA-induced  differentiation  study.  Following  RA  treatment  for  7  days,  Hoxc9+/+ 
spheres  displayed  morphologic  features  of  neuronal  differentiation,  such  as  extensive  neurite  outgrowth  (Fig. 
3A,  Hoxc9+/+).  By  contrast,  most  of  Hoxc9+/'  spheres  showed  no  obvious  morphological  changes  following  RA 
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treatment  (Fig.  3 A,  Hoxc9+/').  At  the  molecular  level,  Hoxc9  heterozygosity  significantly  inhibited  the  ability  of 
RA  to  downregulate  MYCN  and  p75Ngfr,  and  to  upregulate  Gfra3  and  Nefm  (Fig.  3B).  It  has  long  been  known 
that  downregulation  of  MYCN  by  RA  is  critical  for  neuroblastoma  cells  to  differentiate  [19].  p75Ngfr  is  a 
marker  for  neural  crest  stem  cells  [20-22],  which  is  generally  thought  to  give  rise  to  neuroblastoma  stem  cells 
[23,  24],  Together,  these  data  suggest  that  high  Hoxc9  expression  is  essential  for  RA-induced  differentiation  of 
mouse  neuroblastoma  stem  cells. 


Task  3.  To  investigate  the  molecular  mechanism  for  global  control  of  HOXC9-induced  differentiation 
(months  3-32) 

The  goal  of  Task  3  is  to  test  the  hypothesis  that  HOXC9  transcriptionally  activates  the  histone  H3  demethylases 
KDM5B  and  KDM6B  for  global  control  of  its  differentiation  program.  It  is  well  established  that  histone  H3 
methylation  at  specific  lysine  residues  is  a  major  mechanism  for  global  control  of  development  and 
differentiation  programs  by  regulating  the  active  and  silent  states  of  genes  [25-27].  H3K4me3  maintains  genes 
in  an  active  state,  whereas  H3K27me3  keeps  genes  in  a  silent  state.  KDM5B  catalyzes  demethylation  at  H3K4, 
whereas  KDM6B  catalyzes  demethylation  at  H3K27  [28,  29].  We  proposed  that  KDM5B  upregulation  by 
HOXC9  leads  to  a  decrease  in  H3K4me3  levels  at  cell  cycle -promoting  genes,  whereas  KDM6B  upregulation 
by  HOXC9  results  in  a  decrease  in  H3K27me3  levels  at  the  genes  required  for  differentiation. 

The  Task  3  experiments  outlined  in  the  approved  SOW  for  the  second  budget  year  include  KDM6B 
knockdown  studies  and  ChIP-seq  assays. 


Task  3.2. 
a. 


b. 

c. 


KDM5B  and  KDM6B  knockdown  studies  (months  9-20) 


To  examine  the  effects  of  KDM5B  and  KDM6B  knockdown  on  the  induction  of  the  HOXC9- 


induced  differentiation  program,  we  will  infect  BE(2)-C/Tet-Off/HOXC9  cells  cultured  in  the 
presence  of  doxycycline  (Doxy)  with  pLKO.l  lentiviral  constructs  expressing  shRNAs  to  human 
KDM5B  or  KDM6B.  To  examine  the  effects  of  KDM5B  and  KDM6B  knockdown  on  the 


maintenance  of  the  differentiation  program,  BE(2)-C/Tet-Off/HOXC9  cells  will  be  cultured  in  the 
absence  of  Doxy  for  6  days  to  induce  HOXC9  and  differentiation,  followed  by  the  lentiviral 
infection.  pLKO.l-GFP  shRNA  will  be  used  as  negative  control. 

Microarray  assays  will  be  conducted  for  assessing  the  effects  of  KDM5B  and  KDM6B  knockdown 
on  the  ability  of  HOXC9  to  regulate  gene  expression.  There  will  be  a  total  of  18  microarrays:  9 
(GFP-,  KDM5B-  and  KDM6B-shRNA,  3  independent  samples/shRNA)  x  2  (Doxy+/Doxy-)  =18. 
Functional  analyses  will  be  conducted  for  assessing  the  effects  of  KDM5B  and  KDM6B  knockdown 
on  the  ability  of  HOXC9  to  induce  cell  cycle  exit  and  neuronal  differentiation. 


shKDM6B_KDM6  RT-PCR 


Fig.  4.  Compensatory  upregulation 
of  KDM6A  by  KDM6B  knockdown. 
Quantitative  RT-PCR  analysis  of  KDM6A 
and  KDM6B  mRNA  levels  in  BE(2)-C/ 
TetOff/HOXC9  cells  infected  with  lenti- 
viruses  expressing  shRNA  to  GFP  or 
KDM6B. 


As  summarized  in  my  last  annual  report,  we  have  completed  the  study 
to  examine  the  effect  of  KDM5B  knockdown  on  the  ability  of  HOXC9 
to  induce  differentiation. 

During  the  last  budget  year,  we  have  completed  the  study  to  examine 
the  effect  of  KDM6B  knockdown  on  the  ability  of  HOXC9  to  induce 
differentiation.  We  infected  BE(2)-C/Tet-Off/HOXC9  cells  cultured  in 
the  presence  of  doxycycline  (Doxy)  with  lentiviral  constructs 
expressing  shRNA  to  GFP  (control)  or  human  KDM6B.  The  cells 
were  then  cultured  in  the  absence  of  Doxy  for  6  days  to  induce 
HOXC9.  So  far,  we  have  observed  no  significant  effects  of  KDM6B 
knockdown  on  the  ability  of  HOXC9  to  induce  cell  cycle  arrest,  to 
repress  cell  cycle  genes  (i.e.,  CCNA2,  CCNB1  and  CCNE2),  and  to 
induce  neuronal  genes  ( GFRA3  and  NENM)  (data  not  shown). 
However,  we  noticed  that  in  cells  with  KDM6B  knockdown  the 
expression  of  KDM6A  was  upregulated  (Fig.  4).  It  is  known  that 
KDM6A  and  KDM6B  share  the  same  substrate  specificity  [29].  Thus, 


it  is  possible  that  KDM6A  upregulation  masked  the  effect  of  KDM6B  knockdown.  We  have  developed  two 


alternative  strategies  to  address  the  issue. 
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First,  as  proposed  in  my  original  application,  we  have  begun  co-knockdown  experiments  by  sequential 
infection  of  BE(2)-C/Tet-Off/HOXC9  cells  with  lentiviral  constructs  expressing  shRNA  to  KDM6A  and 
KDM6B.  After  verification  by  RT-PCR  and  immunoblotting,  we  will  examine  the  effects  of  KDM6A  and 
KDM6B  co-knockdown  on  the  ability  of  HOXC9  to  induce  differentiation. 

Second,  we  have  taken  the  advantage  of  recent  progress  on  the  development  of  small-molecule  KDM 
inhibitors  and  obtained  a  small-molecule  inhibitor  (GSK-J4)  that  specifically  targets  both  KDM6A  and  KDM6B 
[30].  We  will  treat  BE(2)-C/Tet-Off/HOXC9  cells  with  GSK-J4  and  then  examine  its  effect  on  HOXC9-induced 
differentiation. 

For  microarray  assays,  we  have  already  collected  RNA  samples  from  the  KDM5B  knockdown  study.  We  will 
perform  the  assays  once  we  obtain  RNA  samples  from  the  KDM6  co-knockdown  and  inhibition  studies. 

Task  3.3.  Anti-KDM5B,  -KDM6B,  -H3K4me3,  and  -H3K27me3  ChIP-Seq  assays  (months  21-32): 

a.  ChIP-Seq  assays  will  be  conducted  with  BE(2)-C/Tet-Off/HOXC9  cells  before  (Doxy+)  and  after 
HOXC9  induction  (Doxy-  6  days)  using  ChIP-quality  antibodies  against  KDM5B,  KDM6B, 
H3K4me3,  and  H3K27me3,  followed  by  massively  parallel  DNA  sequencing  (Illumina). 


A 
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Fig.  5.  Epigenetic  regulation  of  transcription  by  HOXC9.  ChIP-seq  tag  profiles  showing  HOXC9  binding  to  representative  target 
genes,  leading  to  changes  in  H3K4me3  and/or  H3K27me3  levels  at  their  loci  (A,  GFRA3;  B,  CCNB1 ;  C,  HOXC  locus). 

We  have  performed  the  proposed  anti-H3K4me3  and  anti-H3K27me3  ChIP-seq  assays.  H3K4me3  maintains 
genes  in  an  active  state,  whereas  H3K27me3  keeps  genes  in  a  silent  state.  Our  ChIP-seq  data  revealed  3  modes 
of  epigenetic  regulation  of  transcription  by  HOXC9  (Fig.  5).  First,  HOXC9-mediated  transcriptional  activation 
of  GFRA3  is  associated  with  a  marked  increase  in  H3K4me3  levels  at  the  GFRA3  promoter  and  a  significant 
reduction  in  H3K27me3  levels  at  the  GFRA3  genomic  locus  (Fig.  5A).  Second,  HOXC9-mediated 
transcriptional  repression  of  CCNB1  (coding  for  cyclin  Bl)  is  associated  with  a  significant  increase  in  H3K27 
levels  at  the  CCNB1  promoter  but  with  no  significant  change  in  H3K4me3  levels  (Fig.  5B).  Third,  HOXC9 
repression  of  anterior  HOXC  genes  is  characterized  by  a  significant  increase  in  H2K27me3  levels  and  a  marked 
reduction  in  H3K4me3  levels  at  their  genomic  loci  (Fig.  5C,  inset).  These  findings  provide  evidence  in  support 
of  our  hypothesis  that  HOXC9  transcriptionally  activates  the  H3K4  demethylase  KDM5B  and  the  H3K27 
demethylase  KDM6B  for  global  control  of  its  differentiation  program. 


3.  KEY  RESEARCH  ACCOMPLISHMENTS: 

•  We  have  completed  the  experiments  proposed  in  MYCN  mouse  tumor  development  studies  and  we  are  in 
the  final  stage  of  data  analysis  to  determine  the  effects  of  Hoxc9  heterozygosity  on  neuroblastoma 
development  and  differentiation. 

•  We  have  developed  a  cell  culture  system  for  enriching  and  long-term  propagating  mouse  neuroblastoma 
sphere-forming  cells  that  possess  cancer  stem  cell  properties. 
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•  We  have  obtained  evidence  that  Hoxc9  heterozygosity  increases  the  population  of  sphere-forming  cells  in 
mouse  neuroblastoma,  suggesting  a  potential  tumor  suppressor  function  of  Hoxc9  in  mouse  neuroblastoma 
development. 

•  We  have  performed  microarray  gene  expression  profding  of  mouse  neuroblastoma  sphere-forming  cells  and 
identified  key  metabolic  pathways  that  are  activated  in  these  cells,  thereby  providing  a  molecular  rational 
for  targeting  these  pathways  for  neuroblastoma  therapy. 

•  We  have  completed  the  study  of  RA- induced  differentiation  of  mouse  neuroblastoma  sphere-forming  cells 
and  obtained  evidence  for  an  important  role  of  Hoxc9  in  promoting  neuroblastoma  stem  cell  differentiation. 

•  We  have  completed  anti-H3K4me3  and  anti-H3K27me3  ChIP-seq  assays  with  BE(2)-C/Tet-Off/HOXC9 
cells  before  (Doxy+)  and  after  HOXC9  induction  (Doxy-  6  days)  and  identified  3  modes  of  epigenetic 
regulation  of  transcription  by  HOXC9. 

4.  CONCLUSION: 

Neuroblastoma  is  a  common  childhood  malignant  tumor  of  the  sympathetic  nervous  system.  Differentiation 
status  in  neuroblastoma  strongly  affects  clinical  outcomes  and  inducing  differentiation  is  a  treatment  strategy  in 
this  disease  [1-5].  We  have  recently  shown  that  HOXC9  is  a  key  regulator  of  neuroblastoma  cell  differentiation 
and  a  prognostic  marker  for  survival  in  neuroblastoma  patients  [6].  Neuroblastoma  differentiation  is 
characterized  at  the  molecular  level  by  repression  of  cell  cycle  genes  and  activation  of  neuronal  differentiation 
markers.  Our  research  during  the  second  budget  year  of  this  award  has  significantly  advanced  our  understanding 
of  the  role  and  the  mechanism  of  HOXC9  action  in  neuroblastoma.  We  developed  a  cell  culture  system  for 
enriching  and  long-term  propagating  mouse  neuroblastoma  sphere-forming  cells  that  possess  cancer  stem  cell 
properties.  With  this  culture  system,  we  were  able  to  show  that  high  Hoxc9  expression  is  essential  for  RA- 
induced  differentiation  of  neuroblastoma  stem  cells.  Moreover,  we  identified  3  modes  of  epigenetic  regulation 
of  transcription  by  HOXC9  that  involve  modulation  of  H3K4me3  and  H3K27me3  levels,  thereby  providing  a 
molecular  basis  for  the  role  of  KDM5B  and  KDM6B  in  the  control  of  neuroblastoma  differentiation.  Finally, 
through  microarray  gene  expression  profiling,  we  obtained  evidence  for  selective  activation  of  the  serine- 
glycine  synthesis  pathway  in  mouse  neuroblastoma  stem  cells.  These  findings  provide  the  first  line  of  direct 
evidence  for  the  presence  of  a  minor  population  of  neuroblastoma  stem  cells  within  primary  tumors  and  suggest 
that  pharmacological  inactivation  of  the  serine-glycine  synthesis  pathway  is  a  promising  therapeutic  strategy  for 
neuroblastoma. 

5.  PUBLICATIONS,  ABSTRACTS,  AND  PRESENTATIONS: 

a)  Manuscripts  and  abstracts 

(1)  Wang,  X.,  Choi,  J.,  Ding,  J.,  Yang,  L.,  Lee,  E.J.,  Ngoka,  Zha,  Y.,  Jin,  B.,  Ren,  M.,  Huang,  S., 
Cowell,  J.,  Shi,  H.,  Cui,  H.  Ding,  H.-F.  HOXC9  directly  regulates  distinct  sets  of  genes  to 
coordinate  diverse  cellular  processes  during  differentiation.  BMC  Genomics  14:830.  2013. 

PMID:  24274069 

(2)  Wang,  X.,  Yang,  L.,  Choi,  J.,  Kitamura,  E.,  Chang,  C.,  Ding,  J.,  Lee,  E.,  Cui,  H.,  Ding,  H.-F. 
Genome-wide  analysis  of  HOXC9-induced  neuronal  differentiation  of  neuroblastoma  cells. 
Genomics  Data  2:50-52,  2014 

(3)  Ding,  J.,  Li,  T.,  Wang,  X.,  Zhao,  E.,  Choi,  J.,  Yang,  L.,  Zha,  Y.,  Dong,  Z.,  Huang,  S.,  Asara, 

J.M.,  Cui,  H.,  Ding,  H.-F.  The  histone  H3  methyltransferase  G9A  epigenetically  activates  the 
serine-glycine  synthesis  pathway  to  sustain  cancer  cell  survival  and  proliferation.  Cell  Metab. 
18:896-907.2013.  PMID:  24315373;  PMCID:  PMC3878056 

b)  Presentations 

(1)  Poster  Presentation:  “Histone  H3  methyltransferase  G9A  epigenetically  activates  the  serine 
synthesis  pathway  to  sustain  cancer  cell  survival  and  proliferation”,  Advances  in 
Neuroblastoma  Research,  Cologne,  Germany.  May  13-16,  2014 
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Abstract 

Background:  Cellular  differentiation  is  characterized  by  the  acquisition  of  specialized  structures  and  functions,  cell 
cycle  exit,  and  global  attenuation  of  the  DNA  damage  response.  It  is  largely  unknown  how  these  diverse  cellular 
events  are  coordinated  at  the  molecular  level  during  differentiation.  We  addressed  this  question  in  a  model  system 
of  neuroblastoma  cell  differentiation  induced  by  HOXC9. 

Results:  We  conducted  a  genome-wide  analysis  of  the  HOXC9-induced  neuronal  differentiation  program.  Microarray 
gene  expression  profiling  revealed  that  HOXC9-induced  differentiation  was  associated  with  transcriptional  regulation  of 
2,370  genes,  characterized  by  global  upregulation  of  neuronal  genes  and  downregulation  of  cell  cycle  and  DNA  repair 
genes.  Remarkably,  genome-wide  mapping  by  ChIP-seq  demonstrated  that  HOXC9  bound  to  40%  of  these  genes, 
including  a  large  number  of  genes  involved  in  neuronal  differentiation,  cell  cycle  progression  and  the  DNA  damage 
response.  Moreover,  we  showed  that  HOXC9  interacted  with  the  transcriptional  repressor  E2F6  and  recruited  it  to  the 
promoters  of  cell  cycle  genes  for  repressing  their  expression. 

Conclusions:  Our  results  demonstrate  that  HOXC9  coordinates  diverse  cellular  processes  associated  with 
differentiation  by  directly  activating  and  repressing  the  transcription  of  distinct  sets  of  genes. 

Keywords:  Neuronal  differentiation,  Cell  cycle  arrest,  DNA  damage  response,  E2F6,  HOXC9,  Neuroblastoma 


Background 

Cellular  differentiation  is  an  essential  process  of  normal 
development  by  which  a  stem  or  progenitor  cell  becomes 
a  post-mitotic,  specialized  cell  with  unique  morphology 
and  function.  In  addition,  it  has  long  been  recognized  that 
differentiated  cells  of  both  normal  and  tumor  origin  are 
defective  in  the  DNA  damage  response  and  repair  at  the 
global  level,  displaying  a  marked  increase  in  sensitivity  to 
ionizing  radiation  and  other  DNA  damaging  agents  [1-3]. 
Consistent  with  these  observations,  recent  studies  have 
shown  that  brain  and  breast  cancer  stem  cells,  a  small 
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subpopulation  of  tumor  cells  thought  to  be  responsible  for 
initiating  and  sustaining  tumor  growth  [4-6],  are  more 
resistant  to  irradiation  and  chemotherapy  than  bulk  tumor 
cells  [7-10].  Particularly  interesting  is  the  observation  that 
inhibition  of  DNA  damage  checkpoint  kinases  can  reverse 
the  radioresistance  of  glioma  stem  cells  [7].  Thus,  a 
molecular  understanding  of  cellular  differentiation  may 
suggest  new  therapeutic  strategies  that  target  both  cell 
proliferation  and  the  DNA  damage  response. 

Among  the  genes  that  have  a  critical  role  in  the  control  of 
cellular  differentiation  are  the  HOX  gene  family  members. 
HOX  genes  encode  a  family  of  transcription  factors  that 
function  as  master  regulators  of  morphogenesis  and  cell  fate 
specification  [11-13].  Dysregulation  of  HOX  gene  expression 
has  been  implicated  in  the  pathogenesis  of  cancers  of  differ¬ 
ent  tissue  types.  In  most  tumor  types,  HOX  genes  function 
as  oncogenes  to  promote  cancer  development  such  as 
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HOXA9  in  leukemia  and  HOXB13  in  ovarian  and  breast 
cancers  [13,14].  However,  in  neuroblastoma,  a  common 
childhood  malignant  tumor  of  the  sympathetic  nervous 
system  [15,16],  there  is  evidence  suggesting  that  HOX  genes 
may  function  as  tumor  suppressors  [13].  Particularly,  down- 
regulation  of  HOXC9  expression  is  significantly  associated 
with  poor  prognosis  in  neuroblastoma  patients  [17,18]. 

Neuroblastoma  cells  can  be  induced  to  undergo  neuronal 
differentiation  by  serum  deprivation  [19],  nerve  growth 
factor  [20]  or  retinoic  acid  (RA)  [21].  RA-induced  neuronal 
differentiation  of  neuroblastoma  cells  is  a  well-established 
model  for  molecular  investigation  of  neuronal  differenti¬ 
ation  [22].  We  recently  reported  that  RA-induced  differenti¬ 
ation  of  neuroblastoma  cells  required  the  activation  of 
several  HOX  genes  [18,23].  Among  them,  HOXC9  appeared 
to  be  a  major  mediator  of  RA  action  in  neuroblastoma  cells. 
HOXC9  expression  was  upregulated  by  RA,  and  silencing 
HOXC9  expression  conferred  resistance  to  RA-induced  dif¬ 
ferentiation.  Importantly,  ectopic  HOXC9  expression  alone 
was  sufficient  to  induce  growth  arrest  and  morphologic  dif¬ 
ferentiation  in  neuroblastoma  cells,  fully  recapitulating  the 
neuronal  differentiation  phenotype  induced  by  RA  [18]. 

Differentiated  neuroblastoma  cells  morphologically  and 
functionally  resemble  mature  peripheral  neurons  charac¬ 
terized  by  G1  arrest,  extensive  neurite  outgrowth,  and  sig¬ 
nificant  resting  potential.  It  has  long  been  observed  that 
differentiated  neuroblastoma  cells  are  highly  sensitive  to 
UV  and  X-ray  radiation  with  a  significantly  reduced  rate 
of  DNA  damage  repair  [20,24-27].  The  molecular  basis  for 
the  differentiation-induced  radiosensitivity  is  not  well 
understood.  The  biological  functions  of  RA  are  mediated 
by  multiple  isotypes  of  RA  receptors  (RARs)  and  retinoid 
X  receptors  (RXRs),  which  form  RAR/RXR  heterodimers 
that  bind  RA  response  elements  in  the  regulatory  regions 
of  RA  target  genes  and  regulate  their  transcription  [28]. 
The  complexity  of  multiple  RARs  and  RXRs  involved  in 
the  action  of  RA  presents  a  daunting  challenge  to  dissect 
the  molecular  mechanism  that  coordinates  the  diverse 
cellular  events  associated  with  differentiation.  Thus, 
the  finding  that  HOXC9  alone  is  able  to  initiate  a  robust 
transcriptional  program  that  drives  neuronal  differentiation 
provides  a  unique  experimental  system  for  this  investiga¬ 
tion.  In  this  study,  we  conducted  genome-wide  profiling  of 
the  HOXC9-initiated  transcriptional  program.  Our  investi¬ 
gation  reveals  that  HOXC9  directly  regulates  the  expression 
of  three  major  sets  of  genes  that  separately  control  neur¬ 
onal  differentiation,  cell  cycle  progression,  and  the  DNA 
damage  response. 

Results 

Gene  expression  profiling  of  HOXC9-induced  neuronal 
differentiation 

To  gain  a  molecular  understanding  of  HOXC9-induced 
differentiation,  we  conducted  microarray  gene  expression 


profiling  of  human  neuroblastoma  BE(2)-C/Tet-Off/myc- 
HOXC9  cells,  which  express  myc-tagged  human  HOXC9 
and  undergo  neuronal  differentiation  in  the  absence  of 
doxycycline  [18]  (Figure  1A).  The  profiling  analysis  identi¬ 
fied  a  total  of  2,370  genes  that  were  differentially  expressed 
(>  +1.5  and  <  -1.5  fold,  P  <0.01),  with  879  genes  being  up- 
regulated  and  1,491  genes  downregulated  (Additional  file  1: 
Table  SI).  Gene  annotation  enrichment  analysis  revealed 
that  HOXC9-induced  differentiation  is  characterized  by  a 
genome-wide  coordination  in  transcriptional  regulation  of 
genes  that  control  neuronal  differentiation,  cell  cycle 
progression,  and  the  DNA  damage  response. 

Global  upregulation  of  neuronal  genes 

Gene  Ontology  (GO)  analysis  of  the  879  HOXC9- 
upregulated  genes  by  DAVID  [29,30]  revealed  that  they 
were  significantly  enriched  for  genes  that  control  nervous 
system  development  such  as  neuron  generation  and  dif¬ 
ferentiation,  axonogenesis,  and  synapse  formation  and 
organization  (Figure  IB  and  Additional  file  2:  Table  S2,  en¬ 
richment  fold  >  2.0,  false  discovery  rate  (FDR)  <1%).  A  total 
of  105  HOXC9-responsive  genes  were  involved  in  nervous 
system  development  (Figure  IB),  accounting  for  approxi¬ 
mately  12%  of  the  879  genes  upregulated  by  HOXC9.  We 
obtained  similar  results  with  Gene  Set  Enrichment  Analysis 
(GSEA),  which  showed  significant  enrichment  of  gene 
sets  involved  in  synaptogenesis  and  neuron  differentiation 
among  the  genes  upregulated  by  HOXC9  (Figure  1C).  Par¬ 
ticularly  significant  was  the  activation  of  ASCII ,  GFRA3 , 
RET,  and  NTN3  (Figure  ID).  ASCF1,  a  member  of  the  basic 
helix-loop-helix  (bHFH)  family  of  transcription  factors,  is  a 
master  regulator  in  the  generation  and  differentiation  of 
sympathetic  neurons  [31,32].  GFRA3  encodes  the  glial  cell 
line-derived  neurotrophic  factor  (GDNF)  family  receptor 
alpha  3  (GFRa3),  which  forms  a  receptor  complex  with 
RET  that  preferentially  binds  the  GDNF  family  ligand  Arte- 
min.  This  receptor  signaling  has  a  critical  role  in  embryonic 
development  of  the  sympathetic  nervous  system,  promoting 
the  survival,  differentiation,  axonal  outgrowth,  and  target 
innervation  of  sympathetic  neurons  [33].  NTN3  (netrin  3) 
belongs  to  a  family  of  extracellular  proteins  that  promote 
axon  growth  and  migration  during  the  development  of  the 
nervous  system  [34].  Ingenuity  Pathways  Analysis  (IPA) 
further  revealed  a  network  of  HOXC9-upregulated  genes 
relevant  to  the  development  and  function  of  sympathetic 
neurons  (Additional  file  3:  Figure  SI).  Together,  these  ana¬ 
lyses  demonstrate  that  HOXC9  activates  a  large  number 
of  neuronal  genes,  providing  the  molecular  mechanism 
for  its  ability  to  induce  neuronal  differentiation  of  neuro¬ 
blastoma  cells. 

Global  downregulation  of  cell  cycle  and  DNA  repair  genes 

GO  analysis  of  the  1,491  HOXC9-downregulated  genes 
revealed  that  they  were  remarkably  enriched  for  genes 
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(See  figure  on  previous  page.) 

Figure  1  Global  upregulation  of  neuronal  genes  by  HOXC9.  (A)  Schematic  of  experiment.  Doxy,  doxycycline.  (B)  DAVID  analysis  of 
upregulated  HOXC9-responsive  genes  for  enriched  GO  biological  process  categories  (enrichment  fold  >  2.0;  FDR  <1%).  The  number  of  genes  for 
each  biological  process  category  is  indicated.  (C)  GSEA  showing  significant  enrichment  of  gene  sets  involved  in  synapse  organization  and  biogenesis 
and  neuron  differentiation  among  the  genes  upregulated  by  HOXC9.  (D)  Heatmap  of  select  neuronal  genes  activated  by  HOXC9. 


that  control  cell  cycle  progression  and  the  DNA  damage 
response  (Figure  2A  and  Additional  file  4:  Table  S3, 
enrichment  fold  >  2.0;  FDR  <1%).  The  analysis  identified 
206  genes  involved  in  cell  cycle  regulation  and  98  genes 
in  the  DNA  damage  response  (Figure  2A).  Similarly, 
GSEA  showed  that  among  the  genes  downregulated  by 
HOXC9,  those  regulating  mitotic  cell  cycle,  DNA  rep¬ 
lication  and  DNA  repair  were  significantly  enriched 
(Figure  2B).  IPA  further  revealed  that  the  downregulated 
genes  include  most  of  cyclin  (CCN)  and  cyclin-dependent 
kinase  (CDI<)  genes,  and  genes  that  control  DNA  rep¬ 
lication,  mitosis,  double-strand  break  (DSB)  repair, 
base  excision  repair  (BER),  nucleotide  excision  repair 
(NER),  mismatch  repair  (MMR),  and  Fanconi  anemia 
(FA)-mediated  repair  (Additional  file  3:  Figures  S2A-S2E 
and  Additional  file  4:  Table  S3).  These  findings  suggest 
that  global  downregulation  of  cell  cycle  and  DNA  repair 
genes  is  the  primary  cause  of  the  cell  cycle  arrest  and 
attenuation  of  the  DNA  damage  response  associated  with 
neuronal  differentiation. 

Genome-wide  mapping  of  HOXC9-binding  sites 

We  next  asked  how  HOXC9  coordinates  the  expression 
of  distinct  sets  of  genes:  the  upregulation  of  genes  critical 
for  nervous  system  development  and  the  downregulation 
of  genes  essential  for  cell  cycle  progression  and  the  DNA 
damage  response.  Mechanistically,  HOXC9  could  function 
through  a  few  master  transcription  factors,  which  in 
turn  regulate  their  own  subsets  of  target  genes  that 
work  together  to  drive  differentiation.  Alternatively, 
HOXC9  could  directly  regulate  distinct  sets  of  genes  to 
coordinate  the  cellular  events  associated  with  differentiation. 
To  test  these  models,  we  conducted  two  independent 
anti-HOXC9  chromatin  immunoprecipitation  assays 
followed  by  massively  parallel  sequencing  of  the  enriched 
DNA  fragments  (ChIP-seq)  for  genome-wide  mapping  of 
HOXC9-binding  sites.  We  identified  a  total  of  29,221 
HOXC9-binding  peaks  with  FDR  less  than  1%  (Figure  3A 
and  Additional  file  5:  Table  S4).  Scatter  plot  analysis 
(Figure  3B,  R  =  0.93,  correlation  coefficient)  and  ChIP-seq 
tag  profiles  (Figure  3C)  demonstrated  that  the  mapping  data 
were  highly  reproducible  between  the  two  independent 
HOXC9  ChIP-seq  samples.  We  next  analyzed  the  distribu¬ 
tion  of  HOXC9-binding  peaks  within  the  genome  that  was 
classified  into  functional  categories  including  promoters 
(within  5  kb  upstream  of  the  transcription  start  site,  TSS), 
S’-untranslated  regions  (S’-UTRs),  exons,  introns,  B’-UTRs, 


downstream  (within  5 -kb  downstream  of  the  gene), 
and  intergenic  regions  (outside  -5  ~  +5  kb  of  genes).  The 
analysis  revealed  that  a  majority  of  HOXC9-binding  peaks 
were  localized  in  introns  (41.2%)  and  intergenic  regions 
(43.4%)  (Figure  3D).  However,  after  normalization  to  the 
size  of  these  functional  regions,  it  became  clear  that 
HOXC9-binding  peaks  were  highly  enriched  in  gene 
promoters  and  5’-UTRs  (Figure  3E  and  F).  Analysis  of 
the  sequences  covered  by  HOXC9-binding  peaks  with  the 
motif-finding  program  MEME  revealed  that  the  most 
enriched  binding  motif  (T/ATTTAT,  E  value  =  1.6e-35) 
corresponds  to  the  Drosophila  Abd-B  motif  (MAO  165.1, 
Figure  3G)  and  is  highly  homologous  to  the  mouse 
Hoxc9-binding  motif  (ATTTAT)  [35].  HOXC9  is  a 
mammalian  ortholog  of  the  Drosophila  Hox  protein 
Abd-B.  Thus,  myc-tagged  HOXC9  binds  to  cognate 
sequences  in  human  neuroblastoma  cells. 

Genome-wide  identification  of  HOXC9  target  genes 

The  ChIP-seq  assay  revealed  that  a  total  of  4,992  genes 
contained  at  least  one  HOXC9-binding  peak  within  5-kb 
upstream  or  downstream  of  their  genomic  loci  (Figure  3A 
and  Additional  file  6:  Table  S5).  We  next  combined  the 
anti-HOXC9  ChIP-seq  data  with  the  HOXC9  microarray 
data  to  generate  a  list  of  genes  that  were  bound  by  HOXC9 
and  whose  expression  levels  were  significantly  changed  as 
a  result  of  HOXC9  induction  (>  +1.5  and  <  -1.5  fold, 
P  <0.01).  The  analysis  revealed  that  954  genes  or  40.3%  of 
the  2,370  HOXC9-responsive  genes  are  direct  targets  of 
HOXC9,  with  445  and  509  genes  being  upregulated  and 
downregulated,  respectively  (Additional  file  7:  Table  S6). 
GO  analysis  of  HOXC9  direct  target  genes  revealed  a 
transcriptional  program  characterized  by  coordinated 
regulation  of  genes  critical  for  neuron  differentiation, 
cell  cycle  progression,  and  the  DNA  damage  response. 

HOXC9  directly  induces  a  large  number  of 
neuronal  genes 

The  only  sets  of  genes  that  were  significantly  enriched 
among  the  upregulated  HOXC9  direct  target  genes  are 
those  exclusively  involved  in  nervous  system  develop¬ 
ment,  particularly  the  generation  and  differentiation  of 
neurons  and  axonogenesis  (Figure  4A  and  Additional  file  8: 
Table  S7,  enrichment  fold  >  2.0,  FDR  <  5%).  The  57 
HOXC9  direct  target  genes  account  for  54.3%  (57/105)  of 
the  HOXC9-responsive  genes  involved  in  nervous  sys¬ 
tem  development  (Figure  IB).  Among  them  are  ASCL1, 
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Figure  2  Global  downregulation  of  cell  cycle  and  DNA  repair  genes  by  HOXC9.  (A)  DAVID  analysis  of  downregulated  HOXC9-responsive 
genes  for  enriched  GO  biological  process  categories  (enrichment  fold  >  2.0;  FDR  <1%).  The  number  of  genes  for  each  biological  process  category 
is  indicated.  (B)  GSEA  showing  significant  enrichment  of  gene  sets  involved  in  mitotic  cell  cycle,  DNA  replication  and  DNA  repair  among  the 
genes  downregulated  by  HOXC9. 


GFRA3,  RET,  and  NTN3.  Figure  4B  shows  the  ChlP- 
seq  tag  profiles  of  HOXC9  binding  to  the  promoter  re¬ 
gions  of  GFRA3,  RET,  and  NTN3.  As  discussed  above, 
these  genes  have  a  critical  role  in  sympathetic  neuro¬ 
genesis  and  axonogenesis. 


HOXC9  directly  represses  a  large  number  of  genes 
essential  for  cell  cycle  progression  and  the  DNA  damage 
response 

GO  analysis  of  the  downregulated  HOXC9  direct  target 
genes  revealed  that  they  were  significantly  enriched  for 
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Figure  3  ChIP-seq  analysis  of  HOXC9  genomic  distribution.  (A)  Summary  of  HOXC9-binding  peaks  and  genes  identified  by  ChIP-seq.  (B-C) 
Scatter  plot  analysis  (B)  and  ChIP-seq  tag  profiles  of  chromosome  21  (C)  showing  a  high  correlation  between  two  independent  HOXC9  ChIP-seq 
samples.  (D)  Pie  chart  showing  genomic  distribution  of  HOXC9-binding  peaks  relative  to  RefSeq  functional  categories  including  promoters  (within 
5  kb  upstream  of  TSS),  5'  UTRs,  exons,  introns,  3'  UTRs,  downstream  (within  5  kb  downstream  of  the  gene),  and  intergenic  regions  (outside  -5  ~  +5  kb 
of  genes).  (E)  Bar  graph  showing  the  distribution  of  HOXC9-binding  peaks  relative  to  the  functional  categories  after  normalization  to  the  size  of  each 
category  in  100  kb.  (F)  Histogram  showing  the  distribution  of  HOXC9-binding  peaks  relative  to  the  nearest  TSS.  Peaks  were  combined  in  100  bp. 

(G)  Web  logo  showing  the  top  enriched  motif  present  in  HOXC9  ChIP-seq  peaks,  which  corresponds  to  the  binding  site  for  Abd-B,  the  Drosophila 
ortholog  of  mammalian  HOX9-13  paralogs. 


genes  that  control  cell  cycle  progression  and  the 
DNA  damage  response  (Figure  5A  and  Additional  file  9: 
Table  S8),  enrichment  fold  >  2.0,  FDR  <  1%).  The  analysis 
identified  52  cell  cycle  genes  that  were  directly  repressed 
by  HOXC9  (Figure  5A),  accounting  for  25.2%  (52/206)  of 
the  HOXC9-responsive  genes  involved  in  cell  cycle 
regulation  (Figure  2A).  It  was  particularly  striking  that 
the  vast  majority  of  the  HOXC9-repressed  cell  cycle 
genes  are  involved  in  the  control  of  the  M  phase  (n  =  25) 
and  DNA  replication  (n  =  21)  (Figure  5A).  Figure  5B  and 
C  show  the  association  of  HOXC9  with  the  promoter 
regions  of  representative  cell  cycle  genes,  including 
CDC4SL  and  MCM3  (DNA  replication),  and  CCNB1  and 
CDCA8  (M  phase).  CDC45L  and  MCM3  are  components 
of  the  replicative  complex  that  catalyzes  DNA  replication 
during  the  S  phase  [36],  while  CDCA8,  also  known  as 
BOREALIN,  is  a  component  of  the  chromosomal  passenger 
complex  essential  for  mitosis  and  cell  division  [37]. 

We  also  identified  32  genes  associated  with  the  DNA 
damage  response  that  were  directly  repressed  by  HOXC9 


(Figure  5A),  accounting  for  32.7%  (32/98)  of  the 
HOXC9-responsive  genes  involved  in  the  DNA  damage 
response.  Figure  5D  shows  the  binding  of  HOXC9  to  the 
promoter  of  FANCM  and  to  both  the  promoter  and 
3’  region  of  FEN1.  FANCM  is  a  component  of  the 
FANCM-FAAP24-MHF  protein  complex  that  binds 
to  DNA  with  interstrand  cross-links  and  is  responsible  for 
recruiting  the  FA  core  complex  to  the  damaged  site  [38]. 
FEN1  (flap  endonuclease  1)  is  essential  for  DNA  replication 
and  repair  by  removing  RNA  and  DNA  5’  flaps  [39]. 

Collectively,  these  findings  suggest  that  HOXC9  directly 
regulates  the  expression  of  distinct  sets  of  genes  to  coord¬ 
inate  the  molecular  and  cellular  processes  characteristic  of 
neuronal  differentiation. 

HOXC9  targets  E2F6  to  the  promoters  of  cell  cycle  genes 

We  next  sought  to  determine  the  molecular  basis  for 
HOXC9  regulation  of  gene  expression  by  identifying 
HOXC9-interacting  proteins.  We  used  a  myc-tag  antibody 
to  isolate  myc-HOXC9  and  its  associated  proteins  from 
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Figure  4  HOXC9  directly  induces  a  large  number  of  neuronal  genes.  (A)  DAVID  analysis  of  upregulated  HOXC9  target  genes  for  enriched 
GO  biological  process  categories  (enrichment  fold  >  2.0;  FDR  <  5%).  The  number  of  genes  for  each  biological  process  category  is  indicated. 

(B)  ChIP-seq  tag  profiles  showing  HOXC9  binding  to  representative  upregulated  HOXC9  target  genes  involved  in  nervous  system  development 
{GFRA3,  top;  RET,  middle;  NTN3,  bottom). 


nuclear  extracts  of  BE(2)-C/Tet'Off/myc-HOXC9  cells  cul¬ 
tured  in  the  absence  of  doxycycline  for  6  days  (Additional 
file  3:  Figure  S3A).  Mass  spectrometric  analysis  of  two 
independent  samples  identified  E2F6  as  a  HOXC9- 


interacting  protein  (Additional  file  3:  Figure  S3B),  a  well 
characterized  transcriptional  repressor  that  plays  a  major 
role  in  repressing  E2F-responsive  genes  essential  for  cell 
proliferation  [40].  It  is  known  that  E2F  family  proteins 
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Figure  5  HOXC9  directly  represses  a  large  number  of  cell  cycle  and  DNA  repair  genes.  (A)  DAVID  analysis  of  down  regulated  H0XC9  target 
genes  for  enriched  GO  biological  process  categories  (enrichment  fold  >  2.0;  FDR  <1%).  The  number  of  genes  for  each  biological  process  category  is 
indicated.  (B-D)  ChIP-seq  tag  profiles  showing  HOXC9  binding  to  representative  downregulated  HOXC9  target  genes  involved  in  DNA  replication 
(B),  mitosis  (C),  and  DNA  repair  (D). 


(E2F1-6)  share  the  same  core  consensus  G/CTTTG/C 
binding  site  [41].  Interestingly,  GSEA  revealed  significant 
enrichment  of  the  E2F-binding  motif  among  the  genes 
downregulated  by  HOXC9  (Additional  file  3:  Figure  S3C). 
Taken  together,  these  observations  suggest  that  E2F6  has 
an  important  role  in  HOXC9-mediated  repression  of  cell 
cycle  genes. 

To  corroborate  the  finding  of  mass  spectrometry,  we 
performed  co-immunoprecipitation  (Co-IP)  experiments 
using  nuclear  extracts  from  BE(2)-C/Tet-Off/myc-HOXC9 
cells  cultured  in  the  absence  of  doxycycline  for  6  days. 
The  myc-tag  antibody,  but  not  control  IgG,  precipitated 
myc-HOXC9  and  E2F6  (Figure  6 A,  left  panel).  Reciprocally, 


an  E2F6  antibody  precipitated  E2F6  and  myc-HOXC9 
(Figure  6 A,  right  panel).  We  next  performed  size-exclusion 
chromatography  using  the  same  nuclear  extracts.  Immuno- 
blot  analysis  revealed  the  presence  of  HOXC9  (~31  kDa)  in 
two  complexes:  the  larger  complex  (peak  at  fraction  20) 
had  an  estimated  molecular  mass  of  -1,800  kDa  and  the 
other  (peak  at  fraction  34)  of  -250  kDa  (Figure  6B).  A 
significant  amount  of  endogenous  E2F6  (-36  kDa)  co-eluted 
with  the  1,800-kDa  HOXC9  complex,  whereas  MEIS2 
(-37-49  kDa),  which  interacts  with  HOX  proteins  and 
functions  as  a  HOX  cofactor  [12],  exclusively  co-eluted 
with  the  250  kDa-HOXC9  complex  (Figure  6B).  Co-IP  ex¬ 
periments  using  pooled  fractions  confirmed  the  association 
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Figure  6  HOXC9  interacts  with  E2F6  and  recruits  it  to  cell  cycle  genes.  (A)  Reciprocal  Co-IP  of  myc-H0XC9  and  E2F6  in  nuclear  extracts 
from  BE(2)-CAet-Off/myc-HOXC9  cells  cultured  in  the  absence  of  doxycycline  (Doxy)  for  6  days.  (B)  Size-exclusion  chromatography  analysis  of 
complexes  containing  myc-HOXC9,  E2F6  or  MEIS2  in  nuclear  extracts  of  BE(2)-CAet-Off/myc-HOXC9  cells  cultured  in  the  absence  of  Doxy  for 
6  days.  (C)  Co-IP  of  myc-HOXC9  and  E2F6  in  pooled  Superose-6  fractions  19-21.  (D)  ChIP-qPCR  analysis  showing  E2F6  binding  to  specific  promoter 
regions  of  the  cell  cycle  genes  CCNB1  and  CDCA8  in  BE(2)-CAet-Off/myc-HOXC9  cells  before  (Doxy+)  and  after  (Doxy-)  HOXC9  induction.  Dashed  lines 
indicate  IgG  control.  Error  bars  represent  SD  (n  =  3).  Data  were  analyzed  with  unpaired,  two-tailed  Student's  f-test  and  p  values  are  indicated. 


of  HOXC9  with  E2F6  within  the  larger  complex 
(Figure  6C). 

To  determine  whether  the  HOXC9-E2F6  interaction 
plays  a  role  in  recruiting  E2F6  to  HOXC9  target  genes 
in  vivo,  we  performed  anti-E2F6  ChIP  using  BE(2)-C/ 
Tet-Off/myc-HOXC9  cells  before  and  after  HOXC9  induc¬ 
tion.  HOXC9  induction  had  no  apparent  effect  on  E2F6 
expression  as  determined  by  microarray  gene  expression 
profiling  (-1.003  fold).  ChIP-qPCR  assay  revealed  that 
E2F6  was  recruited  to  specific  promoter  regions  of  the 
cell  cycle  genes  CCNB1  and  CDCA8  only  after  HOXC9 
induction  (Figure  6D).  By  contrast,  no  significant  binding 
of  E2F6  to  the  NEFM  promoter  was  observed  before  and 


after  HOXC9  induction  (Additional  file  3:  Figure  S4A).  As 
reported  previously,  NEFM  is  a  neuronal  gene  directly 
activated  by  HOXC9  during  differentiation  [18]  (See  also 
Additional  file  3:  Figure  S4B).  Together,  these  data  suggest 
that  elevated  levels  of  HOXC9  facilitate  the  formation  of  a 
repressive  complex  with  E2F6,  which  is  then  recruited  to 
cell  cycle  but  not  neuronal  genes  during  differentiation. 

E2F6  is  essential  for  HOXC9-induced  cell  cycle  arrest  and 
transcriptional  repression  of  cell  cycle  genes 

To  determine  the  functional  significance  of  the  HOXC9- 
E2F6  interaction,  we  examined  the  effect  of  E2F6 
knockdown  on  HOXC9-induced  growth  arrest.  We 
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depleted  E2F6  using  short  hairpin  RNA  (shRNA)  sequences 
targeting  different  coding  regions  of  the  human  E2F6  gene 
(Figure  7A).  Cells  with  E2F6  knockdown  were  highly 
resistant  to  HOXC9-induced  G1  arrest,  showing  continued 
cell  proliferation  (Figure  7B)  and  cell  cycle  progression 
(Figure  7C)  following  HOXC9  induction.  This  was  accom¬ 
panied  by  a  marked  decrease  in  the  population  of  cells  in 
the  G1  phase  and  a  significant  increase  in  the  population 
of  cells  in  the  S  phase  (Figure  7D).  In  addition,  E2F6 
knockdown  largely  abrogated  the  ability  of  HOXC9  to 
repress  cyclin  A2  and  B1  expression,  but  had  no  significant 
effect  on  HOXC9  induction  of  NEFM  (Figure  7E  and  F),  a 
finding  consistent  with  the  observation  of  no  significant 
E2F6  binding  to  the  NEFM  promoter  during  HOXC9- 
induced  differentiation  (Additional  file  3:  Figure  S4A). 
Together,  these  findings  identify  an  essential  and  specific 
role  for  E2F6  in  HOXC9  induction  of  growth  arrest  and 
repression  of  cell  cycle  genes. 

Discussion 

In  this  report,  we  present  evidence  for  a  master  regulator 
of  development  with  the  capacity  to  coordinate  diverse 
cellular  events  characteristic  of  neuronal  differentiation  by 
simultaneously  and  directly  regulating  distinct  sets  of 
genes  (Figure  8).  Through  gene  expression  profiling,  we 
show  that  HOXC9-induced  neuronal  differentiation  is 
characterized  at  the  molecular  level  by  transcriptional 
regulation  of  2,370  genes,  with  global  activation  of  genes 
that  promote  nervous  system  development  and  repression 
of  genes  that  are  essential  for  cell  cycle  progression  and 
the  DNA  damage  response.  Moreover,  through  a  combin¬ 
ation  of  genome-wide  mapping  of  HOXC9  binding  sites 
and  gene  expression  profiling,  we  show  that  HOXC9 
directly  regulates  the  expression  of  954  genes,  -40%  of  the 
2,370  HOXC9-responsive  genes,  including  a  large  number 
of  genes  required  for  neuronal  differentiation,  cell  cycle 
progression  and  the  DNA  damage  response.  Finally,  we 
identify  an  essential  role  for  E2F6  in  HOXC9  repression  of 
cell  cycle  genes  and  induction  of  G1  arrest. 

Our  findings  that  HOXC9  can  both  activate  and 
repress  gene  transcription  are  consistent  with  previous 
observations  from  the  study  of  spinal  cord  development  in 
chick  and  mouse  embryos.  In  the  developing  spinal  cord, 
Hoxc9  functions  as  a  transcription  activator  to  promote  the 
fate  of  preganglionic  motor  column  (PGC)  neurons  [42], 
most  likely  through  its  interaction  with  the  transcription 
factor  FoxPl  [43,44].  However,  Hoxc9  can  also  specify 
the  fate  of  hypaxial  motor  column  (HMC)  neurons  by 
repressing  the  Hox  genes  that  promote  the  switch  of 
HMC  neurons  to  the  lateral  motor  column  (LMC)  neurons 
[35].  Importantly,  our  study  further  demonstrated  that 
within  the  same  population  of  neuroblastoma  cells,  HOXC9 
could  simultaneously  activate  the  genes  that  promote  neur¬ 
onal  differentiation  and  repress  the  genes  that  are  essential 


for  cell  cycle  progression  and  the  DNA  damage  response. 
While  the  molecular  basis  for  the  transcription  activator 
function  of  HOXC9  in  neuroblastoma  cells  remain  to  be 
defined,  we  showed  that  the  ability  of  HOXC9  to  repress 
cell  cycle  genes  depended  on  its  interaction  with  the 
transcription  repressor  E2F6,  a  member  of  the  E2F  family 
of  transcription  factors  that  have  a  critical  role  in  the 
control  of  cell  proliferation  [40]. 

Cellular  differentiation  is  tightly  linked  to  cell  cycle 
exit,  with  the  differentiated  cell  containing  the  G1  content 
of  DNA.  The  molecular  mechanism  that  couples  cell  cycle 
exit  and  differentiation  is  not  well  understood,  although  it 
is  generally  recognized  that  cell  cycle  regulators  influence 
differentiation,  and  cell  fate  determinants  influence  the 
cell  cycle  [45-48] .  A  primary  example  is  the  CDI<  inhibitor 
p27Kipl  as  a  key  regulator  that  links  cell  cycle  exit 
and  differentiation  during  development.  p27Kipl  induces 
G1  arrest  by  associating  with  CDK/cyclin  complexes  and 
inhibits  their  kinase  activity  [49].  Overexpression  of 
p27Xlci,  a  Xenopus  homolog  of  p27Kipl,  in  Xenopus  retina 
glial  progenitor  cells  promotes  both  cell  cycle  exit  and 
differentiation  [50].  Knockout  and  overexpression  studies 
also  demonstrate  an  important  role  of  p27Kipl  in  neuronal 
differentiation  in  the  mouse  cerebral  cortex  by  stabilizing 
Neurogenin  2  [51],  a  proneural  bHLH  transcription  factor 
with  a  central  role  in  cortical  neurogenesis  [52].  On  the 
other  hand,  cell  fate  determinants  can  also  modulate  the 
expression  of  p27Kipl  for  coordinated  regulation  of  cell 
cycle  exit  and  differentiation.  For  instance,  Drosophila 
proneural  bHLH  proteins  cooperate  with  epidermal 
growth  factor  signaling  to  directly  activate  the  transcrip¬ 
tion  of  Dapaco ,  a  homolog  of  p21Cip/p27Kipl,  during  the 
differentiation  of  photoreceptor  cells  [53] . 

Our  findings  suggest  an  alternative  mechanism  for 
coupling  cell  cycle  exit  and  differentiation.  HOXC9  does 
not  regulate  the  expression  of  CDI<  inhibitors,  including 
p27Kipi  and  p21Cip,  and  overexpression  of  either  p27Kipl 
or  p21Cip  fails  to  stop  the  proliferation  of  BE(2)-C  cells 
[18].  Rather,  HOXC9  induces  G1  arrest  by  directly 
repressing  a  large  number  of  genes  essential  for  cell  cycle 
progression  through  the  S  to  M  phases,  including  cyclin 
Bl,  CDCA3,  CDCA8,  BUB  IB,  MCM3  and  MCM8.  This 
transcriptional  repression  function  of  HOXC9  requires 
E2F6.  We  found  that  HOXC9  interacts  with  E2F6  and 
recruits  it  specifically  to  the  promoters  of  cell  cycle  genes. 
E2F6  lacks  a  transactivation  domain  and  functions  as 
a  transcriptional  repressor  for  E2F-responsive  genes 
that  drive  cell  proliferation  [54-58].  Mechanistically,  E2F6 
interacts  with  chromatin  modifiers  with  transcription  re¬ 
pressor  activity  to  establish  a  repressive  chromatin  structure. 
These  chromatin  modifiers  include  the  DNA  methyltrans- 
ferase  Dnmt3b  [59]  and  polycomb-group  (PcG)  proteins 
[60-63].  In  our  study,  we  identified  HOXC9  and  E2F6 
within  a  complex  of  approximately  1,800  kDa.  Whether  this 
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Figure  7  E2F6  is  essential  for  HOXC9  induction  of  G1  arrest  and  repression  of  cell  cycle  genes.  (A)  Immunoblot  analysis  of  E2F6  levels  in 
BE(2)-C^et-Off/myc-HOXC9  cells  infected  with  lentiviruses  expressing  shRNA  against  GFP  or  various  coding  regions  of  E2F6.  E2F6  levels  were 
quantified  against  (3-actin.  (B-D)  Phase  contrast  imaging  and  growth  assay  (B)  and  cell  cycle  analysis  (C,  D)  showing  E2F6  knockdown  abrogated 
HOXC9-induced  growth  arrest.  Error  bars,  SD  (n  =  4).  (E-F)  Immunoblot  analysis  (E)  and  quantification  (F)  showing  that  E2F6  knockdown  abrogated 
HOXC9  repression  of  cyclins,  but  not  HOXC9  induction  of  NEFM.  HOXC9,  CCNA2,  CCNB1  and  NEFM  levels  were  quantified  against  (3-actin  with  the 
protein  levels  in  GFPsh-expressing  cells  cultured  in  the  presence  of  doxycycline  (Doxy+)  were  defined  as  1.0  (dashed  lines).  Error  bars,  SD  (n  =  3). 

Data  in  (D)  and  (F)  were  analyzed  with  unpaired,  two-tailed  Student's  f-test  and  p  values  are  indicated. 
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Figure  8  HOXC9  directly  regulates  the  expression  of  distinct 
sets  of  genes  to  coordinate  diverse  cellular  events  associated 
with  neuronal  differentiation. 


complex  contains  chromatin  modifiers  is  currently  under 
investigation. 

Terminal  cell  differentiation  is  also  tightly  associated  with 
a  global  reduction  in  DNA  damage  repair  activities 
[1-3].  The  underlying  molecular  mechanism  is  not  well 
understood.  It  has  been  reported  that  El  ubiquitin- 
activating  enzyme  can  complement  nucleotide-excision 
repair  deficiency  in  extracts  from  differentiated  macro¬ 
phages,  suggesting  a  role  of  ubiquitination  in  the  con¬ 
trol  of  the  DNA  damage  response  during  differentiation 
[64].  Our  study  revealed  that  in  HOXC9-induced  neur¬ 
onal  differentiation,  attenuation  of  the  DNA  damage  re¬ 
sponse  resulted  from  global  transcriptional  repression 
of  DNA  repair  genes.  This  finding  provides  a  molecular 
mechanism  for  the  long  observed  differentiation-induced 
radiosensitivity  in  neuroblastoma  cells  [20,24-27].  For 
HOXC9-induced  differentiation,  a  total  of  98  genes  with 
functions  in  the  DNA  damage  response  were  significantly 
downregulated.  These  genes  are  involved  in  all  types 
of  DNA  damage  checkpoints  and  repair  pathways.  Import¬ 
antly,  we  show  that  32  of  the  98  genes  are  direct  tar¬ 
gets  of  HOXC9.  Thus,  to  a  large  extent,  HOXC9 
coordinates  neuronal  differentiation  and  attenuation  of 
DNA  repair  activities  by  simultaneously  activating  neuronal 
genes  and  repressing  DNA  repair  genes.  Since  the  DNA 
damage  response  and  DNA  replication  machineries  share 
many  components,  we  speculate  that  the  downregulation 
of  DNA  repair  genes  during  differentiation  is  a  conse¬ 
quence  of  repression  of  cell  cycle  genes,  particularly  those 
involved  in  DNA  replication. 

The  stem  cell  model  of  cancer  attributes  cancer  growth 
to  a  subpopulation  of  cancer  stem  cells.  It  has  been  shown 
recently  that  cancer  stem  cells  are  intrinsically  resistant 
to  ionizing  radiation  and  chemotherapy,  as  a  result  of 
enhanced  checkpoint  activation  and  more  effective  DNA 


damage  repair  [7-10].  Since  differentiation  is  associated  with 
global  downregulation  of  DNA  repair  activities,  a  combin¬ 
ation  of  differentiation-inducing  agents  and  irradiation  or 
chemotherapy  may  prove  to  be  a  more  effective  therapeutic 
strategy  for  targeting  cancer  stem  cells. 

Conclusions 

Using  neuroblastoma  cell  differentiation  as  an  experimen¬ 
tal  system,  we  delineate  a  molecular  mechanism  by  which 
HOXC9  coordinates  diverse  cellular  processes  associated 
with  differentiation  by  directly  activating  and  repressing 
the  transcription  of  distinct  sets  of  genes. 

Methods 

Cell  culture  and  growth  assays 

The  human  neuroblastoma  cell  line  BE(2)-C  (CRL-2268, 
ATCC)  with  Tet-Off  inducible  expression  of  myc-tagged 
human  HOXC9  has  been  described  previously  [18]. 
For  E2F6  knockdown,  BE(2)-C/Tet-Off/myc-HOXC9 
cells  were  infected  with  lentiviruses  expressing  shRNA 
against  E2F6  (TRCN013819,  E2F6sh-2,  TTTCGAGT- 
TA  A  ATA  A  ACC  AG  C;  TRCN013821,  E2F6sh-4,  ATTG 
GTGATGTCATACACTCT;  TRCN018201,  E2F6sh-6,  ATC- 
CAAAGCATCTTCCATTGC;  Thermo  Fisher  Scientific). 
Cells  were  cultured  in  a  1:1  mixture  of  DMEM  and 
Hams  nutrient  mixture  F12  supplemented  with  10% 
fetal  bovine  serum  (Invitrogen-Gibco)  in  the  presence  or 
absence  of  doxycycline.  Cells  were  examined  and  phase 
contrast  images  captured  using  an  Axio  Observer  micro¬ 
scope  and  AxioVision  software  (Carl  Zeiss  Microimaging), 
and  viable  cell  numbers  were  determined  by  trypan 
blue  exclusion  assay.  For  cell  cycle  analysis,  cells  were 
fixed  in  70%  ethanol,  incubated  with  ribonuclease  A 
(Sigma-Aldrich),  and  stained  with  20  pg/ml  propidium 
iodide  (Invitrogen-Gibco).  Samples  were  analyzed  using 
a  FACSCalibur  system  and  ModFitLT  V3.2.1  software 
(BD  Bioscience). 

Microarray  gene  expression  profiling 

Total  RNA  was  isolated  using  Trizol  (Invitrogen)  from  three 
independent  samples  of  BE(2)-C/Tet-Off/myc-HOXC9 
cells  cultured  in  the  presence  or  absence  of  doxycycline 
for  6  days.  RNA  was  measured  and  quality  assessed  by  a 
NanoDrop  spectrophotometer  and  an  Agilent  2100 
Bioanalyzer  (Agilent  Technologies).  Affymetrix  microarray 
analysis  was  performed  using  the  Human  Gene  1.0  ST 
microarray  chip.  Data  were  normalized,  significance 
determined  by  ANOVA,  and  fold  change  calculated 
with  the  Partek  Genomics  Suite  (Partek  Inc.).  Gene 
annotation  enrichment  analysis  was  performed  with 
DAVID  v6.7  [30],  GSEA  [65],  and  IPA  (Ingenuity® 
Systems  www.ingenuity.com)  for  all  significantly  changed 
genes  (>  +1.5  and  <  -1.5  fold,  P  <  0.01). 
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ChIP-seq  and  ChIP-qPCR 

Two  independent  preparations  of  BE(2)-C/Tet-Off/ 
myoHOXC9  cells  cultured  in  the  presence  or  absence 
of  doxycycline  for  6  days  were  used  for  ChIP.  Cross-linked 
chromatin  DNA  was  sheared  through  sonication  and 
immunoprecipitated  using  mouse  anti-myc  tag  (clone  4A6, 
Millipore)  or  mouse  anti-E2F6  (sc-53273,  Santa  Cruz 
Biotechnologies)  according  to  the  published  procedure 
[66].  For  ChIP-seq,  libraries  were  generated  from  ChIP 
genomic  DNA  samples  according  to  the  Illumina 
ChIP-seq  library  construction  procedure,  and  sequenced 
using  Illumina  Genome  Analyzer  IIx  with  a  read  length  of 
36  or  76  bp.  For  ChIP-qPCR,  ChIP  genomic  DNA  samples 
were  assayed  in  triplicate  by  PCR  using  an  iQ5  real-time 
PCR  system  (Bio-Rad)  and  the  following  primer  sets  that 
cover  the  promoter  regions  of  CCNB1  (CCNB1_2  and 
CCNB1_6),  CDCA8  (CDCA8_5P200  and  CDCA8_5P1K), 
and  NEFM  (NEFM_5P1  and  NEFM_5P2):  CCNB1_2: 
CCAGAGAGTTGTTGCAACGAT,  CTGGAGAGCAGT 
GAAGCCAGT;  CCNB1_6:  GGAAGGATTGATCAAACC 
CAG,  AGTCACGGATCCGAAAGAAGG;  CDCA8_5P200: 
GGTATTGCAGAGCCGCCA,  CCTCCCCACCAACCCAC 
C;  CDCA8_5P1I<:  TGGTGCCCATCAGGAGCC,  GGCT 
ATGGGAGTGATAATC;  NEFM_5P1:  GCAGAAAGTAA 
TAAGCAACAA,  CCTGCCTTCTGTAAAGTATTG;  NEFM_ 
5P2:  CCTTTCCTGATTACTTACTGA,  AGGGACTCCAGA 
CCGAAATAG. 

ChIP-seq  data  analysis 

Raw  Illumina  sequencing  reads  from  the  two  independ¬ 
ent  ChIP  replicates  (repl,  GEO  GSM848788  and  rep2, 
GEO  GSM848789)  in  the  FASTQ  format  were  cleaned 
using  in-house  scripts  by  trimming  sequencing  adaptors 
and  low  quality  bases  in  both  ends  (Q  <  67  in  Illumina  1.5). 
Cleaned  sequences  were  then  mapped  to  the  human 
genome  (hgl9)  using  Novoalign  v2.07  for  identifying 
the  reads  that  were  mapped  uniquely  to  a  single  genomic 
locus.  The  identified  reads  from  the  repl  ChIP  sample 
(GEO  GSM848788)  were  used  for  peak  calling  with 
Model-based  Analysis  of  ChIP-Seq  (MACS  vl.4)  [67],  and 
only  those  peaks  with  FDR  <1%  were  compared  with 
RefSeq  genes  in  the  UCSC  genome  browser  and  classified 
into  functional  categories  such  as  promoters,  S’-UTRs, 
exons,  introns,  S’-UTRs,  downstream,  and  intergenic 
regions.  To  measure  the  correlation  of  two  HOXC9 
replicates,  we  used  200  bp  non-overlapping  windows 
where  a  tag  density  is  defined  as  the  number  of  reads 
in  a  window.  We  calculated  Pearson  correlation  coef¬ 
ficient  with  R  >  0.9  being  highly  correlated.  For  motif 
analysis,  we  extracted  100  bp  flanking  sequences  from 
predicted  peak  summits  and  ran  MEME  for  identify¬ 
ing  statistically  overrepresented  motifs.  We  performed 
MAST  to  search  motifs  in  the  peaks  using  the  model 
built  by  MEME. 


Identification  of  HOXC9  target  genes 

Genes  with  HOXC9-binding  peaks  that  are  non-intergenic 
(i.e.,  within  -5  ~  +5  kb  of  genes)  were  defined  as  HOXC9 
target  genes.  To  correlate  HOXC9  binding  to  gene  expres¬ 
sion,  we  combined  the  HOXC9  ChIP-seq  data  with  the 
HOXC9  microarray  data  using  in-house  scripts  to 
generate  a  list  of  the  genes  whose  regulatory  elements 
are  bound  by  HOXC9  and  whose  expression  levels 
are  significantly  changed  (>  +1.5  and  <  -1.5  fold,  P  <  0.01) 
as  the  result  of  HOXC9  induction.  The  significantly 
up-  and  down-regulated  HOXC9  target  genes  were 
then  subjected  to  gene  annotation  enrichment  analysis 
with  DAVID  v6.7,  GSEA,  and  IPA. 


Immunoprecipitation  and  mass  spectrometric  analyses 

BE(2)-C/Tet-Off/myc-HOXC9  cells  were  cultured  in  the 
absence  of  doxycycline  for  6  days  and  nuclear  extracts  were 
prepared  following  the  Dignam  protocol  [68]  except  that 
buffer  C  contained  300  mM  NaCl.  Extracts  from  1  x  107 
cells  were  incubated  with  Protein  A/G  beads  (Invitrogen) 
coated  with  4  pg  mouse  anti-Myc  tag  (clone  4A6,  Millipore) 
or  mouse  IgG  for  overnight  at  4°C.  The  beads  were  washed 
3  times  with  buffer  C  containing  150  mM  NaCl,  dried  in  a 
SpeedVac,  re-suspended  in  a  buffer  containing  8M  urea, 
5  mM  DTT  and  100  mM  ammonium  bicarbonate,  and 
alkylated  with  15  mM  iodoacetamide  for  1  hour.  After 
alkylation,  unreacted  iodoacetamide  was  removed  by 
15  mM  DTT  and  the  urea  concentration  was  diluted 
to  ~1M  with  a  buffer  containing  50  mM  ammonium 
bicarbonate  and  2  mM  CaCl2.  Immunoprecipitated  pro¬ 
teins  were  digested  with  14  ng/pl  sequencing  grade 
trypsin  (Promega)  for  24  hours  at  37°C.  The  digests 
were  desalted  with  a  Micro  Trap  desalting  cartridge 
(Michrom  BioResources),  and  tryptic  peptides  eluted 
with  LC-MS  Solvent  B  (90/10/0.05%:  Acetonitrile/water/ 
heptafluorobutyric  acid)  and  dried  in  a  SpeedVac.  The 
digests  were  analyzed  by  Nano-HPLC  using  a  Nano  Trap 
column  (CL5/61241/00,  Michrom  BioResources)  and  an 
Agilent  1200  Series  Nano  pump  (Agilent  Technologies) 
equipped  with  a  refrigerated  autosampler.  An  Agilent  1200 
Series  Capillary  LC  loading  pump  was  used  to  introduce  the 
sample  onto  a  Captrap  cartridge  for  sample  concentration 
and  de-salting. 

Data-dependent  MS  and  MS/MS  spectra  were  acquired 
on  an  LTQ  Orbitrap  Discovery  (Thermo  Fisher  Scientific) 
using  2  micro-scans,  with  a  maximum  injection  time  of 
200  ms  with  2  Da  peak  isolation  width.  Six  scan  events 
were  recorded  for  each  data  acquisition  cycle.  The  first 
scan  event,  acquired  by  the  FTMS,  was  used  for  full  scan 
MS  acquisition  from  300-2000  m/z.  Data  were  recorded 
in  the  Centroid  mode  only.  The  remaining  five  scan 
events  were  used  for  collisionally  activated  dissociation 
(CAD):  the  five  most  abundant  ions  in  each  peptide  MS 
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were  selected  and  fragmented  to  produce  product-ion 
mass  spectra. 

Database  searching  and  protein  identification 

All  MS/MS  data  were  analyzed  using  Bio  Works  Rev.3.3.1 
SP1  (Thermo  Fisher  Scientific)  and  XITandem  (thegpm. 
org).  SEQUEST  was  set  up  to  search  NCBInr_Homosa- 
piens_05262011.fasta  (221863  entries)  and  the  human.pro- 
tein_RefSeq_01 192012  database  (33376),  and  XITandem 
was  set  up  to  search  subsets  of  the  databases.  SEQUEST 
and  XITandem  were  searched  with  a  fragment  ion  mass 
tolerance  of  0.80  Da  and  a  parent  ion  tolerance  of  10.0 
PPM.  Scaffold  (Proteome  Software)  was  used  to  validate 
MS/MS -based  peptide  and  protein  identifications.  Peptide 
identifications  were  accepted  if  they  could  be  established 
at  greater  than  95.0%  probability  as  specified  by  the  Pep¬ 
tide  Prophet  algorithm  [69].  Protein  identifications  were 
accepted  if  they  could  be  established  at  greater  than  90.0% 
probability  and  contained  at  least  1  identified  peptide. 
Protein  probabilities  were  assigned  by  the  Protein  Prophet 
algorithm  [70].  Proteins  that  contained  similar  peptides 
and  could  not  be  differentiated  based  on  MS/MS  analysis 
alone  were  grouped  to  satisfy  the  principles  of  parsimony. 
Single-peptide  protein  identification  was  accepted  only 
if  the  protein  was  independently  identified  by  both 
SEQUEST  and  XITandem. 

Size-exclusion  chromatography 

Size-exclusion  chromatography  was  performed  with  a 
Superose-6  10/300  GL  column  (24  ml  bed  volume)  and 
an  AKTA  purifier  (GE  Healthcare).  Nuclear  extracts 
(0.5  ml)  were  loaded  onto  the  column  equilibrated  with 
PBS,  and  0.5  ml  fractions  were  collected  and  analyzed. 

Co-immunoprecipitation 

Nuclear  extracts  or  pooled  Sepharose-6  fractions  were 
incubated  with  protein  A/G  beads  coated  with  mouse 
anti-Myc  tag  (clone  4A6),  mouse  anti-E2F6,  or  control 
mouse  IgG  for  2  hours  at  4°C.  After  washing  3  times 
with  PBS,  the  beads  were  suspended  in  standard  SDS 
sample  buffer  and  analyzed  by  immunoblotting. 

Immunoblotting 

Unless  indicated,  all  antibodies  were  from  Santa  Cruz 
Biotechnologies.  Samples  were  suspended  in  SDS  sample 
buffer  and  boiled.  Proteins  were  separated  on  SDS- 
polyacrylamide  gels,  transferred  to  nitrocellulose  mem¬ 
branes,  and  probed  with  the  following  primary  antibodies: 
rabbit  anti-cyclin  A2  (sc-751,  1:200),  rabbit  anti-cyclin  B1 
(sc-752,  1:200),  mouse  anti-myc-tag  (9E10,  hybridoma 
supernatant,  1:10),  rabbit  anti-E2F6  (sc-22823,  1:200), 
mouse  anti-MEIS2  (63-T,  sc-81986,  1:400),  mouse  anti- 
NEFM  (NF-09,  sc-51683,  1:200),  and  rabbit  anti-(3- 
actin  (600-401-886,  Rockland  Immunochemicals,  1:2000). 


Horseradish  peroxidase-conjugated  goat  anti-mouse  and 
goat  anti-rabbit  IgG  were  used  as  secondary  antibodies. 
Proteins  were  visualized  using  a  SuperSignal  West  Pico 
chemiluminescence  kit  (Pierce,  Thermo  Fisher  Scientific) 
and  quantified  with  Image J  (National  Institutes  of  Health). 
For  visualization  and  quantification  with  the  Odyssey 
system,  goat  anti-mouse  IRDye  800,  anti-rabbit  IRDye 
800,  anti-mouse  IRDye  680,  and  anti-rabbit  IRDye  680 
were  used  as  secondary  antibodies  (LI-COR  Biosciences). 

Statistics 

All  quantitative  data  were  analyzed  and  presented  with 
GraphPad  Prism  5.0f  for  Mac  using  unpaired,  two-tailed 
Students  £-test. 
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Induction  of  differentiation  is  a  therapeutic  strategy  in  neuroblastoma,  a  common  pediatric  cancer  of  the  sympa¬ 
thetic  nervous  system.  The  homeobox  protein  HOXC9  is  a  key  regulator  of  neuroblastoma  differentiation.  To  gain 
a  molecular  understanding  of  the  function  of  HOXC9  in  promoting  differentiation  of  neuroblastoma  cells,  we  con¬ 
ducted  a  genome-wide  analysis  of  the  HOXC9-induced  differentiation  program  by  microarray  gene  expression 
profiling  and  chromatin  immunoprecipitation  in  combination  with  massively  parallel  sequencing  (ChIP-seq). 
Here  we  describe  in  detail  the  experimental  system,  methods,  and  quality  control  for  the  generation  of  the  micro¬ 
array  and  ChIP-seq  data  associated  with  our  recent  publication  [1]. 

©  2014  The  Authors.  Published  by  Elsevier  Inc.  This  is  an  open  access  article  under  the  CC  BY-NC-ND  license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 


Specifications 

Organism/cell 

line/tissue 

Homo  sapiens 

Strain(s) 

Human  neuroblastoma  cell  line  BE(2)-C  (CRL-2268,  ATCC) 

Sequencer  or  array  type 

Microarray:  Affymetrix  Human  Gene  1.0  ST 

ChIP-seq:  Illumina  Genome  Analyzer  IIX 

Data  format 

Microarray  raw  data:  CEL  files 

ChIP-seq  raw  data:  FASTQ. files 

Experimental  factors 

BE(2)-C  cells  without  or  with  HOXC9  overexpression 

Experimental  features 

Microarray  gene  expression  profiling  to  identify  genes 
that  are  regulated  by  HOXC9 

ChIP-seq  to  map  genomic  sites  that  are  bound  by 

HOXC9 

Consent 

N/A 

Sample  source  location 

N/A 

Direct  link  to  deposited  data 

Deposited  data  are  available  here:  https://www.ncbi.nlm.nih.gov/ 
geo/queiy/acc.cgi?acc=GSE34422. 
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Experimental  design,  materials  and  methods 

Cell  line 

The  human  neuroblastoma  cell  line  BE(2)-C  was  obtained  from 
ATCC  (CRL-2268,  ATCC).  BE(2)-C  cells  carry  p53  mutation  and  MYCN 
amplification  [2,3],  It  has  been  shown  previously  that  this  cell  line  con¬ 
tains  a  subpopulation  of  neuroblastoma  cells  capable  of  self-renewal 
and  multi-lineage  differentiation  [4,5].  BE(2)-C  cells  can  differentiate 
into  neurons  following  treatment  with  retinoic  acid  (RA)  [4-7]  or  into 
Schwann-like  cells  when  treated  with  bromodeoxyuridine  [4,5,8].  As  a 
result,  BE(2)-C  cells  have  been  used  as  a  model  system  for  investigating 
agents  and  signaling  pathways  that  control  neuroblastoma  cell  differen¬ 
tiation,  which  represents  a  promising  therapeutic  strategy  for  neuro¬ 
blastoma  [9,10]. 

We  have  recently  identified  HOXC9  as  a  downstream  target  gene 
of  RA  and  an  essential  mediator  of  RA  action  in  neuroblastoma 
differentiation:  HOXC9  expression  is  upregulated  by  RA  and  knock¬ 
down  of  HOXC9  expression  confers  resistance  to  RA-induced  differenti¬ 
ation.  Moreover,  HOXC9  induction  can  fully  recapitulate  the  phenotype 
of  RA  treatment  [7].  To  gain  a  molecular  understanding  of  the  mecha¬ 
nism  by  which  HOXC9  induces  the  neuronal  differentiation  of  neuro¬ 
blastoma  cells,  we  generated  BE(2)-C-derived  cells  with  inducible 
expression  of  myc-tagged  human  HOXC9  in  the  absence  of  doxycycline, 
using  the  Retro-X  Tet-Off  Advanced  Inducible  Gene  Expression  System 
(Cion  tech). 


http://dx.doi.org/l  0.101 6/j.gdata.2014.04.002 
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Fig.  1.  Quality  assessment  of  ChIP-seq  data.  (A)  BoxWhisker  type  plot  showing  per-base  sequence  quality  scores  for  read  1  (positive  strand,  left  panel)  and  read  2  (negative  strand,  right 
panel).  Higher  score  indicates  better  quality.  (B)  Quality  assessment  of  sequencing  sufficiency  using  histogram  and  cumulative  fraction  of  base  quality  scores  in  peaks.  (C)  Tag  density  plot 
showing  the  average  of  the  normalized  tag  densities  for  all  RefSeq  genes.  HOXC9-binding  sites  are  highly  enriched  around  TSS.  TTS,  transcription  termination  sites. 


Microarray  and  quality  control 

To  identify  the  genes  that  are  regulated  by  HOXC9,  we  isolated  total 
RNA  from  three  independent  samples  of  BE(2)-C/Tet-Off/myc-HOXC9 
cells  cultured  in  the  presence  or  absence  of  doxycycline  for  6  days 
using  Trizol  (Invitrogen).  The  quantity  and  quality  of  the  RNA  samples 
were  measured  and  assessed  by  a  NanoDrop  spectrophotometer  and 
Agilent  2100  Bioanalyzer  (Agilent Technologies).  Affymetrix  microarray 
analysis  was  performed  using  the  Human  Gene  1.0  ST  microarray  chip. 

The  quality  of  each  CEL  file  was  assessed  using  Affymetrix  Expres¬ 
sion  Console  Software  according  to  the  Affymetrix  standard  protocol 
(Quality  Assessment  of  Exon  and  Gene  1.0  ST  Arrays,  Affymetrix 
White  Paper,  2007).  Relative  log  expression  (RLE)  was  used  to  identify 
outlier  samples.  In  order  to  monitor  labeling  and  hybridization  quality, 
we  used  polyA-control  RNAs  ( Lys ,  Phe,  Thr  and  Dap)  and  bacterial 
spike-in  controls  ( BioB ,  BioC ,  BioD  and  Cre),  respectively. 

CEL  files  were  imported  into  Partek  Genomics  Suit  using  RMA  nor¬ 
malization.  The  probesets  were  annotated  using  the  HuGene-l_0-st- 
vl  Probeset  Annotations  and  Transcript  Cluster  Annotations.  The  differ¬ 
ential  expressions  were  calculated  using  ANOVA  of  the  Partek  package. 


ChIP 

To  identify  the  genomic  sites  that  are  bound  by  HOXC9,  we  per¬ 
formed  two  independent  ChIP-seq  assays  with  BE(2)-C/Tet-Off/myc- 


HOXC9  cells  cultured  in  the  absence  of  doxycycline  for  6  days  according 
to  the  published  procedure  [11].  Briefly,  4  x  1 07  cells  were  fixed  with  1  % 
formaldehyde  for  10  min  and  quenched  with  0.125  M  glycine  for  5  min. 
After  cell  lysis,  cross-linked  chromatin  DNA  was  sheared 
to  approximately  250  bp  by  sonication  (Model  150E  ultrasonic  dis- 
membrator,  Fisher  Scientific),  and  immunoprecipitated  with  10  jug  of 
mouse  anti-myc  tag  (clone  4A6,  Millipore)  and  80  pi  Dynabeads  Protein 
G  (Invitrogen).  The  immunoprecipitated  HOXC9-DNA  complexes  were 
washed  extensively  and  eluted  with  SDS  buffer,  followed  by  incubation 
overnight  at  65  °C  to  reverse  cross-linking.  The  samples  were  then  treat¬ 
ed  sequentially  with  RNase  A  and  proteinase  K  to  degrade  associated 
RNA  and  proteins,  and  ChIP  DNA  was  purified  by  phenol-chloroform 
extraction  and  ethanol  precipitation.  Input  genomic  DNA  was  purified 
from  an  aliquot  of  chromatin  after  sonication.  DNA  concentration  was 
determined  using  a  PicoGreen  dsDNA  quantitation  assay  kit  (Invitrogen). 

ChIP  library  and  quality  control 

ChIP  libraries  were  generated  according  to  the  Illumina  ChIP-seq 
library  construction  procedure.  Briefly,  10  ng  of  ChIP  DNA  was  end 
repaired  with  T4  DNA  polymerase,  Klenow  fragment  and  T4  polynucle¬ 
otide  kinase  (all  from  Enzymatics),  and  an  adenosine  base  was  then 
added  to  the  3'  end  of  the  blunt  phosphorylated  DNA  fragments  by 
Klenow  fragment  (3'  — >  5'  exo-)  (Enzymatics).  This  was  followed  by 
ligation  of  Illumina  genomic  adapters.  Ligated  DNA  around  the  size  of 
250-300  bp  was  isolated  by  electrophoresis  through  a  3%  NuSieve  3:1 


52 


X.  Wang  et  al.  /  Genomics  Data  2  (2014)  50-52 


CL 


Top  10  terms  for  cellular  metabolic  processes_HOXC9  downregulated  genes 


purine  base  biosynthetic  process 
deoxyribonucleotide  biosynthetic  process 
nucleobase  biosynthetic  process 
serine  family  amino  acid  biosynthetic  process 
ribonucleoside  monophosphate  metabolic  process 
cholesterol  biosynthetic  process 
sterol  biosynthetic  process 
cellular  amino  acid  biosynthetic  process 
glutamine  family  amino  acid  metabolic  process 
amine  biosynthetic  process 


0 


2 


4 


6 


8 


Fold  enrichment  (p  <  0.001) 


Fig.  2.  HOXC9-induced  neuroblastoma  cell  differentiation  is  characterized  by  reduced  production  of  biosynthetic  precursors.  Gene  ontology  analysis  of  downregulated  HOXC9-responsive 
genes  reveals  the  enrichment  of  genes  involved  in  the  biosynthesis  of  nucleotides,  sterol,  and  amino  acids  (enrichment  fold  >2.0,  p  <  0.001 ). 


agarose  gel  and  amplified  by  PCR  using  Phusion  DNA  Polymerase  (Ther¬ 
mo  Scientific).  The  PCR  products  were  then  purified  using  an  Agencourt 
Ampure  kit  (Beckman  Coulter)  to  remove  primer  dimmers. 

Agilent  2100  Bioanalyzer  was  employed  to  examine  the  normal 
size  distribution  and  possible  primer  or  linker  contamination  of  each 
library.  The  library  concentration  was  measured  using  both  Qubit 
dsDNA  BR  Assay  Kit  (Invitrogen)  and  real-time  qPCR  method.  Real¬ 
time  qPCR  was  performed  on  StepOne  Plus  Real-Time  PCR  systems 
(Applied  Biosystems)  using  KAPA  SYBR  FAST  universal  qPCR  kit 
(KAPA  Biosystems).  qPCR  primers,  1.1  and  2.1,  were  designed  based 
on  Illumina  Sequencing  Library  qPCR  Quantification  Guide.  PhiX  Con¬ 
trol  v3  (Illumina)  was  used  in  the  reaction  to  plot  the  standard  curve. 

ChIP-seq  data  quality  control 

ChIP  DNA  was  sequenced  using  Illumina  Genome  Analyzer  IIX  with  a 
read  length  of  36  (paired-end).  We  performed  several  analyses  to  assess 
the  quality  of  the  sequencing  data.  Box  whisker  type  plot  analysis  shows 
per  base  sequence  quality  scores  for  read  1  (Fig.  1  A,  left  panel)  and  read 
2  (Fig.  1A,  right  panel).  Red  and  blue  lines  represent  the  median  and 
mean  value,  respectively;  the  yellow  box  indicates  the  inter-quartile 
range  (25-75%);  and  the  upper  and  lower  whiskers  indicate  the  10% 
and  90%  points  (FASTQC)  [12].  Fig.  IB  shows  the  quality  assessment  of 
sequencing  sufficiency  using  histogram  and  cumulative  fraction  of 
base  quality  scores  in  peaks.  We  also  analyzed  tag  density  to  visualize 
the  quality  of  sequencing  and  mapping,  and  to  examine  the  global 
profile  of  HOXC9  ChIP-seq  based  on  genes,  which  shows  that  HOXC9- 
binding  peaks  are  highly  enriched  around  transcription  start  sites 
(TSS)  (Fig.  1C). 

Discussion 

We  described  here  a  dataset  composed  of  microarray  gene  ex¬ 
pression  profiling  of  HOXC9-responsive  genes  and  ChIP-seq  data  for 
genome-wide  mapping  of  HOXC9-binding  sites  in  the  human  neuro¬ 
blastoma  BE(2)-C  cells.  With  this  dataset,  we  were  able  to  demonstrate 
that  HOXC9  directly  controls  and  coordinates  the  expression  of  a  large 
number  of  genes  involved  in  neuron  genesis  and  differentiation,  cell 
cycle  progression,  and  DNA  damage  repair  [1  ].  We  believe  that  this 
dataset  would  be  particularly  valuable  for  investigating  the  cellular 
processes  associated  with  differentiation  and  underlying  molecular 


mechanisms.  For  example,  we  found  that  HOXC9-induced  differentia¬ 
tion  is  associated  with  a  marked  reduction  in  cellular  metabolism,  as  a 
result  of  global  downregulation  of  genes  involved  in  the  biosynthesis 
of  nucleotides,  sterol,  and  amino  acids  (Fig.  2).  This  finding  sheds 
new  light  on  the  interplay  between  metabolic  reprogramming  and 
cancer  cell  differentiation,  and  suggests  that  reprogramming  cellular 
metabolism  may  represent  a  promising  strategy  for  promoting  cancer 
cell  differentiation. 


Conflict  of  interest 

The  authors  have  no  conflicts  of  interest. 


References 

[1]  X.  Wang,  J.H.  Choi,  J.  Ding,  L.  Yang,  L.C.  Ngoka,  E.J.  Lee,  Y.  Zha,  L.  Mao,  B.  Jin,  M.  Ren,  J. 
Cowell,  S.  Huang,  H.  Shi,  H.  Cui,  H.F.  Ding,  HOXC9  directly  regulates  distinct  sets  of 
genes  to  coordinate  diverse  cellular  processes  during  neuronal  differentiation. 
BMC  Genomics  14  (2013)  830. 

[2]  D.  Tweddle,  A.  Malcolm,  N.  Bown,  A.  Pearson,  J.  Lunec,  Evidence  for  the  development 
of  p53  mutations  after  cytotoxic  therapy  in  a  neuroblastoma  cell  line.  Cancer  Res.  61 
(2001)  8-13. 

[3]  N.  Keshelava,  J.J.  Zuo,  P.  Chen,  S.N.  Waidyaratne,  M.C.  Luna,  C.J.  Gomer,  T.J.  Triche,  C. 
P.  Reynolds,  Loss  of  p53  function  confers  high-level  multidrug  resistance  in  neuro¬ 
blastoma  cell  lines.  Cancer  Res.  61  (2001)  6185-6193. 

[4]  H.  Cui,  J.  Ma,  J.  Ding,  T.  Li,  G.  Alam,  H.F.  Ding,  Bmi-1  regulates  the  differentiation  and 
clonogenic  self-renewal  of  I-type  neuroblastoma  cells  in  a  concentration-dependent 
manner.  J.  Biol.  Chem.  281  (2006)  34696-34704. 

[5]  RA.  Ross,  B.A.  Spengler,  C.  Domenech,  M.  Porubcin,  W.J.  Rettig,  J.L.  Biedler,  Human 
neuroblastoma  I-type  cells  are  malignant  neural  crest  stem  cells.  Cell  Growth  Differ. 
6  (1995)  449-456. 

[6]  N.  Sidell,  Retinoic  acid-induced  growth  inhibition  and  morphologic  differentiation 
of  human  neuroblastoma  cells  in  vitro.  J.  Natl.  Cancer  Inst.  68  (1982)  589-596. 

[7]  L.  Mao, J.  Ding,  Y.  Zha,  L.  Yang,  BA.  McCarthy,  W.  King,  H.  Cui,  H.F.  Ding,  HOXC9  links 
cell-cycle  exit  and  neuronal  differentiation  and  is  a  prognostic  marker  in  neuroblas¬ 
toma.  Cancer  Res.  71  (2011)  4314-4324. 

[8]  T.  Sugimoto,  T.  Kato,  T.  Sawada,  Y.  Horii, J.T.  Kemshead,  T.  Hino,  H.  Morioka,  H.  Hosoi, 
Schwannian  cell  differentiation  of  human  neuroblastoma  cell  lines  in  vitro  induced 
by  bromodeoxyuridine.  Cancer  Res.  48  (1988)  2531-2537. 

[9]  C.P.  Reynolds,  K.K.  Matthay,  J.G.  Villablanca,  B.J.  Maurer,  Retinoid  therapy  of  high-risk 
neuroblastoma.  Cancer  Lett.  197  (2003)  185-192. 

[10]  S.L.  Volchenboum,  S.L.  Cohn,  Progress  in  defining  and  treating  high-risk  neuroblastoma : 
lessons  from  the  bench  and  bedside.  J.  Clin.  Oncol.  27  (2009)  1003-1004. 

[11]  T.I.  Lee,  S.E.  Johnstone,  R.A.  Young,  Chromatin  immunoprecipitation  and  microarray- 
based  analysis  of  protein  location.  Nat.  Protoc.  1  (2006)  729-748. 

[12]  S.  Andrews,  FastQC:  a  quality  control  tool  for  high  throughput  sequence  data,  http:// 
www.bioinformatics.babraham.ac.uk/projects/fastqc/  2011. 


Cel 


Cell  Metabolism 

Article 


The  Histone  H3  Methyltransferase  G9A  Epigenetically 
Activates  the  Serine-Glycine  Synthesis  Pathway 
to  Sustain  Cancer  Cell  Survival  and  Proliferation 


Jane  Ding,1  -1 0  Tai  Li,1  -6-1 0  Xiangwei  Wang,7  Erhu  Zhao,1  >6  Jeong-Hyeon  Choi,1  >2  Liqun  Yang,6  Yunhong  Zha,  Zheng  Dong,3 
Shuang  Huang,4  John  M.  Asara,9  Hongjuan  Cui,6  and  Han-Fei  Ding1’4-5* 

1  Cancer  Center 

department  of  Biostatistics  and  Epidemiology 
department  of  Cellular  Biology  and  Anatomy 
department  of  Biochemistry  and  Molecular  Biology 
department  of  Pathology 

Medical  College  of  Georgia,  Georgia  Regents  University,  Augusta,  GA  30912,  USA 

6State  Key  Laboratory  of  Silkworm  Genome  Biology,  Institute  of  Sericulture  and  System  Biology,  Southwest  University, 

Chongqing  400716,  China 

department  of  Urology,  Second  Affiliated  Hospital,  Third  Military  Medical  University,  Chongqing  400037,  China 
department  of  Neurology,  The  First  Hospital  of  Yichang,  Three  Gorges  University  College  of  Medicine,  Yichang  423000,  China 
division  of  Signal  Transduction,  Beth  Israel  Deaconess  Medical  Center,  Harvard  Medical  School,  Boston,  MA  02115,  USA 
10These  authors  contributed  equally  to  this  work 
Correspondence:  hding@gru.edu 
http://dx.doi.Org/1 0.1 01 6/j.cmet.201 3.1 1 .004 


SUMMARY 

Increased  activation  of  the  serine-glycine  biosyn¬ 
thetic  pathway  is  an  integral  part  of  cancer  meta¬ 
bolism  that  drives  macromolecule  synthesis  needed 
for  cell  proliferation.  Whether  this  pathway  is  under 
epigenetic  control  is  unknown.  Here  we  show  that 
the  histone  H3  lysine  9  (H3K9)  methyltransferase 
G9A  is  required  for  maintaining  the  pathway  enzyme 
genes  in  an  active  state  marked  by  H3K9  monome- 
thylation  and  for  the  transcriptional  activation  of 
this  pathway  in  response  to  serine  deprivation.  G9A 
inactivation  depletes  serine  and  its  downstream 
metabolites,  triggering  cell  death  with  autophagy  in 
cancer  cell  lines  of  different  tissue  origins.  Higher 
G9A  expression,  which  is  observed  in  various 
cancers  and  is  associated  with  greater  mortality  in 
cancer  patients,  increases  serine  production  and 
enhances  the  proliferation  and  tumorigenicity  of 
cancer  cells.  These  findings  identify  a  G9A-depen- 
dent  epigenetic  program  in  the  control  of  cancer 
metabolism,  providing  a  rationale  for  G9A  inhibition 
as  a  therapeutic  strategy  for  cancer. 


INTRODUCTION 

Histone  lysine  methylation  has  a  central  role  in  the  control  of 
gene  transcription.  The  histone  lysine  methylation  state  is 
controlled  by  histone  lysine  methyltransferases  (KMTs)  and 
demethylases  (KDMs).  The  enzymatic  reactions  catalyzed  by 
KMTs  and  KDMs  depend  on  metabolic  coenzymes  including 
S-adenosylmethionine  (SAM),  flavin  adenine  dinucleotide 


(FAD),  and  a-ketoglutarate  (a-KG).  KMTs  catalyze  lysine  methyl¬ 
ation  using  SAM  as  the  methyl  group  donor,  whereas  LSD 
(KDM1A  and  KDM1B)  and  JmjC  domain-containing  KDMs 
(KDM2-KDM8)  require  FAD  and  a-KG  for  demethylation,  respec¬ 
tively  (Black  et  al.,  2012;  Mosammaparast  and  Shi,  2010).  Their 
dependence  on  metabolic  coenzymes  suggests  that  KMTs  and 
KDMs  could  reprogram  gene  expression  in  response  to  changes 
in  cellular  metabolism.  This  notion  has  also  led  to  the  provocative 
hypothesis  that  KMTs  and  KDMs  may  contribute  to  metabolic 
control  through  transcriptional  regulation  (Teperino  et  al.,  2010). 

G9A,  also  known  as  EHMT2,  is  a  H3K9  methyltransferase 
that  has  a  primary  role  in  catalyzing  monomethylation  and  dime- 
thylation  of  H3K9  (H3K9me1  and  H3K9me2)  in  euchromatin 
(Peters  et  al.,  2003;  Rice  et  al.,  2003;  Shinkai  and  Tachibana, 
2011;  Tachibana  et  al.,  2002),  with  H3K9me1  being  associated 
with  transcriptional  activation  and  H3K9me2  with  transcriptional 
repression  (Black  et  al.,  2012;  Mosammaparast  and  Shi,  2010). 
Elevated  levels  of  G9A  expression  have  been  observed  in 
many  types  of  human  cancers,  and  G9A  knockdown  has  been 
shown  to  inhibit  the  proliferation  of  cancer  cell  lines  (Chen 
et  al.,  2010;  Cho  et  al.,  2011;  Huang  et  al.,  2010;  Kondo  et  al., 
2008).  The  molecular  basis  of  G9A  action  in  the  control  of  cancer 
cell  proliferation  is  not  well  understood.  In  this  study,  we  identify 
an  essential  role  of  G9A  in  sustaining  cancer  cell  survival  and 
proliferation  by  transcriptional  activation  of  the  serine-glycine 
biosynthetic  pathway.  Our  findings  provide  direct  evidence  for 
a  G9A-dependent  epigenetic  program  in  the  control  of  amino 
acid  production  and  cancer  metabolism. 

RESULTS 

G9A  Is  Essential  for  Sustaining  Cancer  Cell  Proliferation 
and  Survival 

We  examined  the  role  of  G9A  in  cell  survival  and  proliferation  in 
human  cancer  cell  lines  of  different  tissue  origins,  including  the 
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Figure  1.  G9A  Inhibition  or  Silencing 
Induces  Autophagy 

(A)  Electron  micrographs  of  U20S  cells  with  or 
without  5  fiM  BIX  for  24  hr. 

(B)  Immunofluorescence  of  LC3B-positive  auto¬ 
phagosomes  in  U20S  cells  with  or  without  5  [xM 
BIX  for  24  hr. 

(C)  Immunoblotting  of  LC3B  in  U20S  cells  with  or 
without  5  jxM  BIX  or  10  [xM  chloroquine  (CQ)  for 
24  hr.  LC3B-II  levels  were  quantified  against 
a-tubulin. 

(D)  Immunoblotting  of  G9A  and  LC3B  in  U20S  and 
SHEP1  cells  expressing  GFP  shRNA  (GFP-sh)  or 
individual  G9A-sh  sequences.  G9A  and  LC3B-II 
levels  were  quantified  against  (3-actin. 

(E)  Fluorescence  and  quantification  of  LC3B- 
positive  autophagosomes  in  U20S/mCherry- 
EGFP-LC3B  cells  with  or  without  G9A  knockdown. 
Error  bars  represent  SD  (n  =  1 00  cells).  **p  <  0.001 . 
See  also  Figures  SI  and  S2. 
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bladder  (J82),  bone  (U20S),  brain  (U251),  breast  (MCF10A  and 
MCF7),  cervix  (HeLa),  colon  (FICT116  and  RKO),  liver  (Hep2G), 
lung  (HI 299),  and  sympathetic  nervous  system  (BE(2)-C,  SMS- 
KCNR,  and  SHEP1).  We  treated  these  cell  lines  with  BIX01294 
(BIX),  a  small  molecule  inhibitor  of  G9A  (IC50  =  1 .7  fiM)  (Kubicek 
et  al.,  2007).  BIX  at  2-5  pM  significantly  reduced  the  global 
levels  of  FI3K9me1  and  H3K9me2  (Figure  SI  A  available  online) 
and  completely  inhibited  the  proliferation  of  all  the  cancer 
cell  lines  examined  (Figure  SIB  for  representative  cell  lines). 
In  addition,  we  observed  a  significant  decrease  in  cell  survival 
following  BIX  treatment  (Figure  SIC,  BIX-5  pM_5d).  To  confirm 
that  BIX  targets  G9A  to  inhibit  cell  proliferation  and  survival,  we 
examined  the  effect  of  G9A  silencing  by  small  hairpin  RNA 
(shRNA).  G9A  knockdown  exerted  a  pronounced  inhibitory 
effect  on  cell  proliferation  and  survival  (Figures  S1D-S1G). 
Together,  these  findings  indicate  an  essential  role  of  G9A  in 
sustaining  cell  proliferation  and  survival  in  a  wide  range  of  cancer 
cell  lines. 

G9A  Inhibition  or  Silencing  Induces  Autophagy 

An  early  and  prominent  morphological  feature  of  the  cells  with 
G9A  inhibition  or  silencing  was  the  appearance  of  numerous 
cytoplasmic  vesicles  and  vacuoles  (Figures  SIC  and  S1G) 
that  morphologically  resemble  autophagosomes,  a  double- 
membraned  structure  that  sequesters  cellular  organelles, 
proteins,  and/or  lipids  during  autophagy.  Thus,  we  examined 
the  possibility  that  G9A  inhibition  or  silencing  might  induce 


autophagy  by  electron  microscopy  for 
ultrastructural  morphology,  by  immuno¬ 
blotting  for  detecting  the  lipidation  of 
LC3  (microtubule-associated  protein  light 
chain  3),  and  by  immunofluorescence  for 
monitoring  the  formation  of  LC3-positive 
puncta.  LC3  is  the  mammalian  homolog 
of  the  yeast  autophagy-related  protein 
Atg8  and  is  proteolytically  processed  to 
LC3-I  by  the  Atg4  protease  following 
translation.  Upon  autophagy  induction, 
LC3-I  is  converted  to  the  lipidated  LC3-II  form,  which  is  then 
incorporated  into  the  autophagosomal  membrane,  resulting  in 
the  redistribution  of  LC3  from  a  diffuse  pattern  to  a  punctate 
pattern.  Mammalian  cells  express  three  LC3  isoforms  (LC3A, 
LC3B,  and  LC3C),  with  LC3B-II  levels  correlating  with  the 
steady-state  levels  of  autophagosomes  (Klionsky  et  al.,  2008; 
Mizushima  et  al.,  2010). 

Electron  microscopy  revealed  numerous  double-membraned 
vacuoles  in  BIX-treated  cells  that  contained  fragments  of 
the  endoplasmic  reticulum  and  other  cytoplasmic  components 
(Figure  1A).  Immunofluorescence  confirmed  that  these  vesicles 
and  vacuoles  were  LC3B  positive  (Figure  IB).  To  determine 
whether  the  accumulation  of  autophagosomes  resulted  from 
an  increase  in  autophagosome  formation  or  from  a  block  in 
autophagosome  turnover,  we  examined  BIX-induced  LC3B-II 
production  in  the  presence  or  absence  of  chloroquine  (CQ), 
which  blocks  lysosomal  acidification  and  thus  degradation 
of  autophagosomal  components.  The  presence  of  CQ  further 
increased  LC3B-II  levels  in  BIX-treated  cells  (Figure  1C), 
indicating  that  BIX  increased  LC3B-II  production  rather  than 
blocking  its  degradation. 

Dose-titration  and  time  course  studies  revealed  that  different 
tumor  cell  lines  varied  in  their  sensitivity  to  BIX.  For  the  most 
sensitive  SHEP1  and  U20S  cells,  BIX  as  low  as  1  (iM  triggered 
significant  levels  of  autophagy,  as  determined  by  LC3B-II  pro¬ 
duction  (Figure  S2A).  The  maximal  level  of  autophagy  occurred 
at  6  hr  (SHEP1  and  U2QS)  or  24  hr  (HeLa)  following  5  jiM  BIX 
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treatment  (Figure  S2B).  We  further  confirmed  the  ability  of  BIX  to 
induce  autophagy  with  immortalized  Atg5+/+  and  Atg5~'~  mouse 
embryonic  fibroblasts  (MEFs).  Atg5  is  required  for  autophago¬ 
some  formation  (Kuma  et  al.,  2004;  Mizushima  et  al.,  2001). 
BIX  at  the  concentrations  of  2-5  fiM  was  able  to  induce  LC3B- 
positive  puncta  in  Atg5+/+,  but  not  in  Atg5~^,  MEFs  (Figure  S2C). 

To  verify  that  BIX  targets  G9A  to  induce  autophagy,  we  exam¬ 
ined  the  effect  of  BIX  on  cells  with  G9A  overexpression.  Both 
SHEP1  and  U20S  cells  overexpressing  G9A  were  highly  resis¬ 
tant  to  BIX-induced  autophagy,  as  evidenced  by  a  significant 
reduction  in  LC3B-II  levels  following  BIX  treatment  compared 
to  control  cells  (Figure  S2D).  Additionally,  cells  with  G9A 
silencing  fully  recapitulated  the  autophagy  phenotype  induced 
by  BIX,  showing  a  significant  increase  in  the  production  of 
LC3B-II  (Figure  ID)  and  LC3B-positive  puncta  (Figure  IE). 
Collectively,  these  results  demonstrate  that  G9A  inhibition  or 
silencing  induces  cell  death  with  autophagy. 

G9A  Inhibition  or  Silencing  Transcriptionally  Represses 
Serine-Glycine  Biosynthesis 

To  investigate  the  molecular  basis  of  autophagy  induction  by 
G9A  inhibition,  we  examined  the  kinase  activity  of  the  mamma¬ 
lian  target  of  rapamycin  complex  1  (mTORCI),  which  has  a  major 
role  in  the  control  of  autophagy  by  integrating  a  wide  range  of 
signals  including  concentrations  of  nutrients  and  growth  factors 
(Fie  and  Klionsky,  2009;  Kim  and  Guan,  201 1 ;  Rabinowitz  and 
White,  2010).  mTORCI  inhibits  autophagy  by  suppressing  the 
autophagy-activating  kinase  ULK1.  Nutrient  or  growth  factor 
deprivation  inhibits  mTORCI  activity,  thereby  activating  ULK1 
and  inducing  autophagy.  Immunoblotting  revealed  that  BIX 
significantly  inhibited  mTORCI  kinase  activity  as  demonstrated 
by  a  marked  decrease  in  the  phosphorylation  of  the  ribosome 
protein  S6  kinase  (S6K)  (Figure  S3A),  a  key  mTORCI  substrate. 
This  finding  suggests  that  G9A  inhibition  triggers  autophagy 
by  interfering  with  cell  growth  or  survival  signals  upstream  of 
mTORCI . 

Given  the  role  of  G9A  in  epigenetic  control  of  transcription,  we 
performed  microarray  profiling  to  identify  potential  G9A  targets 
involved  in  autophagy  induction.  A  total  of  615  BIX-responsive 
genes  (>±1 .50-fold,  p  <  0.01)  were  identified,  with  302  genes 
being  upregulated  and  313  downregulated  (Table  SI).  Gene 
Ontology  (GO)  analysis  revealed  that  among  the  genes  down- 
regulated  by  BIX,  those  within  the  serine-glycine  biosynthetic 
pathway  were  significantly  enriched  (Figure  2A),  including  phos- 
phoglycerate  dehydrogenase  (PFIGDFI),  phosphoserine  amino¬ 
transferase  1  (PSAT1),  phosphoserine  phosphatase  (PSPH), 
and  serine  hydroxymethyltransferase  2  (SHMT2)  (Figure  2B). 
We  confirmed  the  ability  of  BIX  to  downregulate  mRNA  expres¬ 
sion  of  these  genes  by  quantitative  reverse-transcription  PCR 
(qRT-PCR)  in  four  different  cancer  cell  lines  (Figures  2C  and 
S3B).  Time  course  studies  revealed  significant  downregulation 
of  these  genes  within  2  and  6  hr  of  BIX  treatment  in  U20S  and 
HeLa  cells,  respectively  (Figure  2C),  which  occurred  at  least 
several  hours  before  most  of  the  cells  underwent  autophagy 
(Figure  S2B).  We  noticed  that  SHMT1,  though  downregulated 
in  HeLa,  SHEP1 ,  and  U20S  cells,  was  activated  in  BE(2)-C  cells 
following  BIX  treatment  (Figure  S3B).  The  significance  of  this 
apparently  cell-type-dependent  regulation  of  SHMT1  expression 
is  currently  under  investigation.  In  agreement  with  the  microarray 


and  qRT-PCR  data,  we  observed  a  significant  decrease  in 
PHGDH  protein  levels  in  BIX-treated  cells  (Figure  S3C).  Similarly, 
G9A  silencing  resulted  in  marked  downregulation  of  the  same 
group  of  genes  at  both  mRNA  and  protein  levels  (Figures  2D, 
S3D,  and  S3E). 

Chromatin  immunoprecipitation  and  quantitative  PCR  (ChlP- 
qPCR)  assay  revealed  that  BIX  treatment  significantly  reduced 
the  H3K9me1  levels  around  the  transcriptional  start  sites 
(TSSs)  of  PHGDH  and  PSAT1  (Figure  2E,  H3K9me1).  By 
contrast,  H3K9me2  levels  in  the  same  regions  were  either 
unchanged  or  significantly  increased  following  BIX  treatment 
(Figure  2E,  H3K9me2).  Together,  these  results  indicate  that  the 
serine-glycine  biosynthetic  pathway  is  under  the  direct  tran¬ 
scriptional  control  of  G9A,  primarily  through  the  regulation  of 
H3K9me1  levels  associated  with  the  pathway  enzyme  genes. 

Supplemental  Serine  Rescues  the  Cell  Death  Phenotype 
of  G9A  Inhibition 

The  serine-glycine  biosynthetic  pathway  is  an  important  source 
of  precursors  and  coenzymes  for  the  biosynthesis  of  amino 
acids,  nucleotides,  and  lipids  (de  Koning  et  al.,  2003;  Kalhan 
and  Hanson,  201 2)  and  is  crucial  for  cancer  cell  survival  and  pro¬ 
liferation  (Locasale,  2013;  Locasale  et  al.,  2011;  Pollari  et  al., 
2011;  Possemato  et  al.,  2011).  In  addition,  amino  acid  depriva¬ 
tion  is  a  major  trigger  of  autophagy  (Kim  and  Guan,  2011). 
Thus,  we  investigated  the  possibility  that  suppression  of  this 
biosynthetic  pathway  might  be  a  cause  of  the  cell  death  pheno¬ 
type  induced  by  G9A  inhibition  or  silencing.  Gas  chromatog¬ 
raphy-mass  spectrometry  (GC-MS)  analysis  revealed  that  only 
serine  and  glycine  levels  were  significantly  reduced  in  U20S 
cells  within  4  hr  of  BIX  treatment  (Figure  3A),  demonstrating 
that  G9A  activity  is  essential  for  maintaining  the  intracellular 
steady-state  levels  of  serine  and  glycine.  We  further  assessed 
the  activity  of  this  biosynthetic  pathway  by  [U-13C]  glucose  flux 
analysis  using  liquid  chromatography-tandem  mass  spectrom¬ 
etry  (LC-MS/MS).  BIX  treatment  significantly  decreased  the 
incorporation  of  [U-13C]  glucose  into  3-phosphoserine  and 
serine  (Figure  3B).  Importantly,  addition  of  serine  to  the  culture 
medium  significantly  diminished  the  effect  of  BIX  (Figures  3C 
and  3D)  or  G9A  silencing  (Figure  3E)  on  cell  proliferation  and 
autophagy  in  all  the  cancer  cell  lines  examined.  By  contrast, 
supplemental  glycine  had  only  a  small  protective  effect  on 
BIX-treated  cells,  and  addition  of  both  serine  and  glycine  was 
no  more  effective  than  adding  serine  alone  (Figures  3F,  S4A, 
and  S4B).  We  confirmed  these  findings  with  cell-permeable 
methyl-serine-ester  and  methyl-glycine-ester  (Figure  3G).  Other 
individual  amino  acid  supplements  all  failed  to  prevent  cell  death 
with  autophagy  induced  by  BIX  (Figures  3F,  S4A,  and  S4C), 
providing  further  evidence  for  the  specificity  of  serine  action. 

The  observation  that  serine,  but  not  glycine,  was  able  to 
rescue  the  cell  death  phenotype  suggests  that  the  production 
of  5,10-methylenetetrahydrofolate  (5,10-MTHF)  might  be  crucial 
for  cell  survival  and  proliferation.  The  interconversion  between 
serine  and  glycine  is  catalyzed  by  SHMT:  serine  conversion  to 
glycine  generates  5,10-MTHF,  whereas  glycine  conversion  to 
serine  consumes  5,10-MTHF  (Figure  3H).  To  test  this  idea,  we 
generated  SHEP1  cells  overexpressing  SHMT1,  a  cytoplasmic 
enzyme,  or  SHMT2,  a  predominantly  mitochondrial  enzyme 
(Figure  31).  Their  overexpression  alone  had  no  significant  effect 


898  Cell  Metabolism  18,  896-907,  December  3,  2013  ©2013  Elsevier  Inc. 


Cell  Metabolism 

G9A  Control  of  Serine  Metabolism 


Cell 


Downregulated  BIX-responsive  genes_GO  terms 


L-serine  biosynthesis 
L-serine  metabolism 
serine  family  amino  acid  biosynthesis 
serine  family  amino  acid  metabolism 
cellular  response  to  unfolded  protein 
neutral  amino  acid  transport 
cellular  amino  acid  biosynthesis 
amino  acid  activation 
cell  cycle  arrest 
RTK  signaling  pathway  Jm 
regulation  of  RNA  metabolic  process  . 
receptor  protein  signaling  pathway  . 

regulation  of  transcription _ 


B  3-Phosphoglycerate 

||  PHGDH 

3-Phosphohydroxypyruvate 

It? 


PSAT1 
.9 

3-Phosphoserine 
PSPH 


I  p,! 


.8 


Serine 


it? 


SHMT2 
8 


10 


20  30  40 

Fold  enrichment 


50 


60 


Glycine 


Un 

BIX_2h 

BIX_6h 

BIX_24h 

BIX_48h 


U20S_BIX 


HeLa  BIX 


Un 

BIX_6h 

BIX_24h 


U20S_G9Ash  HeLa_G9Ash 

■ 

1.5n 


H3K9me2 


■  GFP-sh 
□  G9A-sh4 


<3^  XV  <X  an  XV 


,,  V 


Figure  2.  G9A  Inhibition  or  Silencing  Represses  Serine-Glycine  Biosynthesis 

(A)  GO  analysis  of  downregulated  BIX-responsive  genes  (enrichment  fold  >2.0,  p  <  0.01). 

(B)  Serine-glycine  biosynthetic  pathway  with  the  indicated  fold  changes  in  mRNA  expression  determined  by  microarray. 

(C  and  D)  qRT-PCR  analysis  of  mRNA  expression  of  the  pathway  enzyme  genes  in  U20S  and  HeLa  cells  with  or  without  5  fiM  BIX  (C)  or  G9A  silencing  (D). 
Error  bars  represent  SD  (n  =  3). 

(E)  ChIP-qPCR  analysis  of  H3K9me1  and  H3K9me2  levels  at  the  promoters  of  PHGDH  and  PSAT1  in  U20S  cells  with  or  without  5  BIX  for  4  hr.  Data  on  the 
GAPDH  coding  region  are  shown  as  positive  control.  Dashed  line  indicates  IgG  control.  Error  bars  represent  SD  of  triplicate  and  are  representatives  of  two 
independent  experiments. 

*p  <  0.01 ,  **p  <  0.001 .  See  also  Figure  S3,  Table  SI ,  and  Table  S3. 


on  BIX-induced  cell  death  (Figure  3J).  However,  overexpression 
of  SHMT2,  but  not  SHMT1 ,  synergized  with  supplemental  serine 
to  enhance  cell  survival  and  proliferation  in  the  presence  of  BIX 
(Figure  3J).  This  observation  is  consistent  with  our  model,  since 
previous  studies  have  suggested  that  SHMT2  has  a  major  role  in 


the  conversion  of  serine  to  glycine  and  5,10-MTHF,  whereas 
SHMT1  primarily  catalyzes  the  conversion  of  glycine  to  serine 
(Herbig  et  ah,  2002;  Narkewicz  et  ah,  1996;  Pfendner  and  Pizer, 
1980;  Stover  et  al.,  1997;  Tibbetts  and  Appling,  2010;  Yoshida 
and  Kikuchi,  1970).  Collectively,  these  results  indicate  that 
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Figure  3.  Supplemental  Serine  Rescues  the 
Phenotype  of  G9A  Inhibition  or  Silencing 

(A)  GC-MS  analysis  of  amino  acid  levels  in  U20S 
cells  with  or  without  5  jxM  BIX  for  4  hr.  Error  bars 
represent  SD  (n  =  6). 

(B)  LC-MS/MS  analysis  of  [U-13C]  glucose  flux  into 
3-phosphoserine  (pSer)  and  serine  (Ser)  biosyn¬ 
thesis  in  U20S  cells  with  or  without  5  BIX  for 
1 6  hr.  The  fraction  of  labeled  to  total  pSer  and  Ser, 
13C/(13C  +  12C),  was  calculated  and  normalized  to 
control.  Error  bars  represent  SD  (n  =  3). 

(C)  Growth  assay  of  SHEP1  and  U20S  cells  with  or 
without  5  [iM  BIX  or  supplemental  serine  (Ser). 
Error  bars  represent  SD  (n  =  4). 

(D)  Immunoblotting  of  LC3B  in  the  cell  samples 
from  (C).  LC3B-II  levels  were  quantified  against 
a-tubulin. 

(E)  Growth  assay  of  the  indicated  cell  lines  with 
G9A  silencing  in  the  presence  or  absence  of  sup¬ 
plemental  serine  for  4  days.  Dashed  line  indicates 
the  number  of  cells  plated  at  time  zero.  Error  bars 
represent  SD  (n  =  4). 

(F  and  G)  Growth  assay  of  U20S  cells  in  the 
presence  of  5  [xM  BIX  with  or  without  the  indi¬ 
cated  supplemental  amino  acids  for  2  days.  Gly, 
glycine;  Gly-ester,  methyl-glycine-ester;  Ser, 
serine;  Ser-ester,  methyl-serine-ester.  Dashed  line 

(F)  indicates  the  number  of  cells  plated  at  time 
zero.  Error  bars  represent  SD  (n  =  4). 

(FI)  Schematic  of  the  interconversion  between 
serine  and  glycine.  THF,  tetrahydrofolate;  5,10- 
MTHF,  5,1 0-methylenetetrahydrofolate. 
(l)qRT-PCR  analysis  of  SHMT  mRNA  expression  in 
SHEP1  cells  overexpressing  SHMT1  or  SHMT2. 
Error  bars  represent  SD  (n  =  3). 

(J)  Growth  assay  of  SHEP1  cells  expressing  GFP, 
SHMT1 ,  or  SHMT2  in  the  presence  of  5  ^M  BIX 
with  or  without  supplemental  Gly-ester  or  Ser- 
ester.  Dashed  line  indicates  the  number  of  cells 
plated  at  time  zero.  Error  bars  represent  SD  (n  =  3). 
*p  <  0.01 ,  **p  <  0.001 .  See  also  Figure  S4. 
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maintaining  the  production  of  serine  and  its  downstream  metab¬ 
olites,  including  5,10-MTHF,  is  a  major  mechanism  by  which 
G9A  sustains  cancer  cell  survival  and  proliferation. 

G9A  Is  a  Key  Component  of  the  Serine  Deprivation 
Response,  Linking  Serine  Production,  Ribosome 
Biogenesis,  and  Cell-Cycle  Progression 

It  has  been  shown  recently  that  serine  deprivation  leads  to  tran¬ 
scriptional  activation  of  the  serine  biosynthetic  pathway  (Ye 
et  al.,  2012;  Figure  4A).  In  light  of  our  findings  above,  we  asked 
whether  G9A  has  a  role  in  activation  of  this  pathway  in  response 
to  serine  deprivation.  Addition  of  BIX  completely  abrogated  the 
induction  of  PHGDH  and  PSAT1  by  serine  deprivation  (Figures 
4B  and  4C),  supporting  a  physiological  function  of  G9A  in  medi¬ 
ating  the  serine  deprivation  response. 

Amino  acid  availability  controls  the  rates  of  protein  synthesis 
and  cell  growth  (Kilberg  et  al.,  2005),  and  it  has  been  reported 


recently  that  G9A  is  required  for  activation 
of  rRNA  transcription  (Yuan  et  al.,  2007),  a 
rate-limiting  step  in  ribosome  biogenesis. 
We  therefore  asked  whether  G9A  has  a  role  in  linking  serine 
metabolism  to  ribosome  biogenesis  and  cell  growth.  To  test 
this  idea,  we  performed  temporal  transcriptome  profiling  of 
BIX-responsive  genes  by  RNA  sequencing  (RNA-seq),  which 
revealed  distinctive  patterns  of  gene  expression  in  cells  that 
were  treated  with  BIX  for  6  or  24  hr  (Figure  4D).  In  agreement 
with  the  microarray  data  shown  above  (Figure  2A;  Table  SI), 
GO  analysis  of  the  genes  downregulated  by  BIX  at  6  hr  revealed 
significant  enrichment  of  genes  within  the  serine  biosynthetic 
pathway  (Figure  4E;  Table  S2).  However,  cells  treated  with  BIX 
for  24  hr  displayed  a  gene  expression  pattern  characterized  by 
significant  downregulation  of  genes  that  control  cell-cycle  pro¬ 
gression  (Figure  4F),  including  cyclins  A2  and  B2,  and  CDC25C 
(Table  S2).  We  confirmed  the  RNA-seq  result  by  qRT-PCR, 
which  showed  sequential  downregulation  of  serine  synthesis 
and  cell-cycle  genes  following  BIX  treatment  (Figure  4G).  We 
also  confirmed  the  time-dependent  downregulation  of  cyclins 
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Figure  4.  G9A  Links  Serine  Production, 
Ribosome  Biogenesis,  and  Cell  Proliferation 

(A)  qRT-PCR  analysis  of  PHGDH  and  PSAT1  mRNA 
levels  in  SHEP1  cells  cultured  in  serine-deficient 
(-Ser)  medium  for  0,  2,  6,  and  24  hr.  Error  bars 
represent  SD  (n  =  3). 

(B)  qRT-PCR  analysis  of  PHGDH  mRNA  levels  in 
SHEP1  cells  cultured  in  -Ser  medium  with  or 
without  5  pM  BIX  or  supplemental  serine  for  0,  2,  6, 
and  24  hr.  Error  bars  represent  SD  (n  =  3). 

(C)  Immunoblotting  of  PHGDH  and  PSAT1  in 
SHEP1  cells  cultured  in  control  (with  serine) 
or  -Ser  medium  with  or  without  5  jxM  BIX  or  sup¬ 
plemental  serine  for  24  hr.  PHGDH  and  PSAT1 
levels  were  quantified  against  a-tubulin. 

(D)  Heatmap  showing  differential  gene  expression 
in  U20S  cells  without  (control)  or  with  5  pM  BIX 
for  6  (BIX6h)  or  24  hr  (BIX24h). 

(E  and  F)  GO  analysis  of  downregulated  BIX- 
responsive  genes  at  the  6  hr  (E)  or  24  hr  (F)  time 
point.  Shown  are  top  ten  GO  biological  process 
terms. 

(G  and  H)  qRT-PCR  (G)  and  immunoblot  (H) 
analyses  of  cell-cycle  genes  in  SHEP1  cells  with  or 
without  5  pM  BIX  for  the  indicated  time.  Error  bars 
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(I)  GO  analysis  of  the  genes  downregulated  by 
BIX  at  24  hr  versus  6  hr.  Shown  are  top  ten  GO 
biological  process  terms. 

See  also  Table  S2  and  Table  S3. 
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A2  and  B1  by  immunoblotting  (Figure  4H).  Importantly,  we  found 
marked  downregulation  of  genes  involved  in  ribosome  biogen¬ 
esis  in  cells  treated  with  BIX  for  24  hr  compared  to  those  treated 
for  6  hr  (Figure  41;  Table  S2),  a  finding  consistent  with  the 
reported  role  of  G9A  in  activation  of  rRNA  transcription  (Yuan 
et  al.,  2007).  Collectively,  these  results  suggest  a  key  role  of 
G9A  in  coupling  serine  sensing  to  the  transcription  of  genes 
that  control  protein  synthesis  and  cell  proliferation. 

G9A  Promotes  Cell  Proliferation  and  Tumorigenicity 

G9A  is  overexpressed  in  various  human  cancers  (Chen  et  al., 
2010;  Cho  et  al.,  2011;  Huang  et  al.,  2010;  Kondo  et  al.,  2008), 


and  it  has  been  reported  recently  that 
elevated  levels  of  G9A  expression  posi¬ 
tively  correlate  with  disease  progression 
and  poor  prognosis  in  lung  cancer  (Chen 
et  al.,  2010).  In  line  with  these  reports, 
we  found  that  higher  G9A  expression 
was  significantly  associated  with  reduced 
overall  survival  in  neuroblastoma  patients 
(Figure  5A).  These  observations  prompted 
us  to  investigate  the  functional  conse¬ 
quence  of  high  G9A  expression.  Ectopic 
expression  of  G9A  in  SHEP1  and  U20S 
cells  significantly  increased  cell  prolifera¬ 
tion  (Figure  5B)  and  the  expression  of 
genes  that  promote  cell-cycle  progres¬ 
sion  (Figure  5C),  indicating  that  high  G9A 
expression  alone  is  sufficient  to  confer 
a  growth  advantage  to  cancer  cells.  More¬ 
over,  G9A  overexpression  markedly  enhanced  the  anchorage- 
independent  growth  (Figure  5D)  and  tumorigenicity  (Figures  5E 
and  5F)  of  SHEP1  and  U20S  cells,  demonstrating  that  G9A 
has  transforming  potential.  Together,  our  results  from  both 
human  primary  tumors  and  cell  lines  suggest  an  oncogenic 
function  of  G9A  in  cancer  development. 

Activation  of  the  Serine-Glycine  Biosynthetic  Pathway 
Is  Essential  for  the  Oncogenic  Activity  of  G9A 

Given  the  critical  role  of  the  serine-glycine  biosynthetic  pathway 
in  sustaining  cancer  cell  survival  and  proliferation  (Locasale 
et  al.,  201 1 ;  Pollari  et  al.,  201 1 ;  Possemato  et  al.,  2011),  we  asked 
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Figure  5.  G9A  Promotes  Cell  Proliferation 
and  Tumorigenicity 

(A)  Kaplan-Meier  overall  survival  for  G9A/EHMT2 
expression  in  primary  neuroblastoma  tumors 
(Versteeg  NB88  data  set). 

(B)  Growth  assay  of  SHEP1  and  U20S  cells  with  or 
without  G9A  overexpression.  Error  bars  represent 
SD  (n  =  4). 

(C)  qRT-PCR  analysis  of  cyclin  A2  and  B1  mRNA 
levels  in  SHEP1  cells  with  or  without  G9A  over¬ 
expression.  Error  bars  represent  SD  (n  =  3).  **p  < 
0.001. 

(D)  Soft  agar  assay  of  SHEP1  and  U20S  cells  with 
or  without  G9A  overexpression.  Shown  are  rep¬ 
resentatives  of  three  independent  experiments. 

(E  and  F)  Xenograft  assay  of  SHEP1  and  U20S 
cells  with  or  without  G9A  overexpression.  Images 

(E)  were  taken  on  day  42  after  inoculation.  Tumor 
weight  (F)  was  analyzed  by  scatterplot  with  hori¬ 
zontal  lines  indicating  the  mean. 


whether  elevated  levels  of  G9A  lead  to  increased  activation  of 
this  pathway  and  whether  this  is  a  major  mechanism  for  the 
oncogenic  activity  of  G9A.  Ectopic  expression  of  G9A  in 
SHEP1  and  U20S  cells  markedly  increased  the  mRNA  (Fig¬ 
ure  6A)  and  protein  (Figure  6B)  levels  of  the  pathway  enzyme 
genes.  Moreover,  ChIP  in  combination  with  DNA  sequencing 
(ChIP-seq)  revealed  a  significant  increase  in  H3K9me1  levels 
at  these  gene  loci  (Figure  6C  for  PHGDH  and  SHMT2).  We 
observed  no  significant  change  in  H3K9me2  levels  in  the  same 
regions  (Figure  6C).  These  results  suggest  that  G9A  activates 
the  expression  of  the  pathway  enzyme  genes  by  increasing  the 
H3K9me1  levels  associated  with  their  loci.  As  expected,  G9A 
overexpression  significantly  increased  the  intracellular  levels  of 
serine  and  glycine  (Figure  6D),  and  the  flux  of  glucose  into  the 
serine  biosynthetic  pathway  (Figure  6E). 

To  determine  whether  the  oncogenic  function  of  G9A  depends 
on  its  ability  to  activate  the  serine-glycine  pathway,  we  silenced 
the  expression  of  PHGDH  or  PSAT1  in  G9A-overexpressing 
U20S  cells  by  shRNA  (Figures  7A  and  7B).  Knockdown  of  either 
PHGDH  or  PSAT1  abolished  the  ability  of  G9A  to  promote  cell 
proliferation  in  a  manner  dependent  on  levels  of  PHGDH  or 
PSAT1  downregulation  (Figures  7C  and  7D).  In  fact,  G9A-over- 
expressing  cells  with  more  than  90%  of  PHGDH  or  PSAT1 
knockdown  failed  to  proliferate  in  the  absence  of  supplemental 
serine  (Figures  7C  and  7D).  We  obtained  essentially  the  same 
results  with  G9A-overexpressing  SHEP1  (Figures  S5A-S5D). 
Collectively,  these  data  demonstrate  that  G9A  promotes  cell 
proliferation  by  transcriptionally  activating  the  serine-glycine 
biosynthetic  pathway. 

The  serine  biosynthetic  genes  are  direct  transcriptional  targets 
of  activating  transcription  factor  4  (ATF4)  (Adams,  2007;  Seo 
et  al.,  2009;  Ye  et  al.,  2012),  which  is  activated  in  response  to 
amino  acid  deprivation  (Kilberg  et  al.,  2009).  It  has  been  shown 
that  ATF4  knockdown  or  deficiency  led  to  lower  expression  of 
serine  biosynthetic  enzymes  (Ye  et  al.,  2012),  suggesting  an 
essential  role  of  ATF4  in  maintaining  transcriptional  activation 
of  this  pathway.  To  examine  the  role  of  ATF4  in  G9A  action,  we 


silenced  ATF4  expression  in  G9A-overexpressing  U20S  (Fig¬ 
ure  7E)  or  SHEP1  cells  (Figures  S5E),  which  resulted  in  marked 
downregulation  of  PHGDH  and  PSAT1  at  both  mRNA  and  pro¬ 
tein  levels  (Figures  7E  and  7F  and  S5E  and  S5F).  ATF4  silencing 
also  significantly  reduced  the  proliferation  of  G9A-overexpress- 
ing  U20S  (Figure  7G)  or  SHEP1  cells  (Figure  S5G).  Together, 
these  results  indicate  that  ATF4  is  required  for  G9A  to  transcrip¬ 
tionally  activate  serine  biosynthesis  and  to  stimulate  cell 
proliferation. 

DISCUSSION 

In  this  study,  we  show  that  the  histone  H3K9  methyltransferase 
G9A  is  essential  for  transcriptional  activation  of  the  serine- 
glycine  biosynthetic  pathway  by  specifically  marking  the 
pathway  enzyme  genes  with  H3K9me1,  an  epigenetic  marker 
associated  with  active  chromatin  (Black  et  al.,  2012;  Mosamma- 
parast  and  Shi,  201 0).  We  further  present  evidence  that  G9A  is  a 
component  of  the  molecular  pathway  that  couples  serine 
sensing  to  the  transcriptional  control  of  serine  production,  ribo¬ 
some  biogenesis,  and  cell  proliferation.  Finally,  we  provide 
evidence  that  G9A  has  an  oncogenic  function  in  cancer  develop¬ 
ment  by  conferring  survival  and  growth  advantages  to  tumors 
through  increasing  the  production  of  serine  and  its  downstream 
metabolites.  Our  study  uncovers  an  epigenetic  mechanism  for 
the  control  of  amino  acid  metabolism  and  provides  a  molecular 
explanation  for  the  functional  significance  of  G9A  overexpres¬ 
sion  observed  in  various  human  cancers  (Chen  et  al.,  2010; 
Cho  et  al.,  201 1 ;  Huang  et  al.,  201 0;  Kondo  et  al.,  2008). 

Cell  growth  and  proliferation  depend  on  protein  synthesis, 
which  must  in  turn  match  the  cellular  nutritional  status  including 
amino  acid  availability  (Grummt  and  Ladurner,  2008;  Kilberg 
et  al.,  2005).  The  mechanistic  connection  between  these  cellular 
processes  remains  to  be  clarified.  Our  study  suggests  an  epige¬ 
netic  mechanism  for  the  control  of  ribosome  biogenesis  and  cell 
proliferation  in  response  to  serine  abundance.  We  show  that 
G9A  is  essential  for  the  serine  deprivation  response  that  leads 
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Figure  6.  G9A  Transcriptionally  Activates 
Serine-Glycine  Biosynthesis 

(A)  qRT-PCR  analysis  of  mRNA  levels  of  serine- 
glycine  biosynthetic  enzyme  genes  in  SHEP1  and 
U20S  cells  with  or  without  G9A  overexpression. 
Error  bars  represent  SD  (n  =  3). 

(B)  Immunoblotting  of  PHGDH  and  PSAT1  in 
SHEP1  and  U20S  cells  with  or  without  G9A 
overexpression.  PHGDH  and  PSAT1  levels  were 
quantified  against  a-tubulin. 

(C)  ChIP-seq  tag  profiles  for  H3K9me1  and 
H3K9me2  levels  in  the  PHGDH  and  SHMT2  loci  in 
U20S  cells  with  or  without  G9A  overexpression. 
Horizontal  bars  represent  chromatin  regions  in 
which  H3K9  methylation  levels  are  significantly 
changed  (G9A/GFP),  and  the  numbers  indicate 
fold  enrichment. 

(D)  GC-MS  analysis  of  serine  and  glycine  levels  in 
U20S  cells  with  or  without  G9A  overexpression. 
Error  bars  represent  SD  (n  =  6). 

(E)  LC-MS/MS  analysis  of  [U-13C]  glucose  flux 
into  3-phosphoserine  (pSer)  and  serine  (Ser) 
biosynthesis  in  U20S  cells  with  or  without  G9A 
overexpression.  The  fraction  of  labeled  to  total 
pSer  and  Ser,  13C/(13C  +  12C),  was  calculated 
and  then  normalized  to  GFP  control.  Error  bars 
represent  SD  (n  =  3). 

*p  <  0.05,  **p  <  0.001 .  See  also  Table  S3. 


to  transcriptional  activation  of  genes  for  serine  biosynthesis. 
Moreover,  we  show  that  G9A  activity  or  expression  levels  affect 
serine  production  and  cell  proliferation.  Specifically,  G9A  inhibi¬ 
tion  initiated  a  transcriptional  program  that  results  in  sequential 
repression  of  genes  required  for  serine  production,  ribosome 
biogenesis,  and  cell-cycle  progression,  leading  to  lower  levels 
of  serine  biosynthesis  and  cell  growth  arrest.  Conversely,  G9A 
overexpression  activated  genes  for  serine  biosynthesis  and 
cell-cycle  progression,  leading  to  higher  levels  of  serine  produc¬ 
tion  and  cell  proliferation.  These  results  are  of  particular  interest 
in  the  context  of  recent  work  showing  that  G9A  is  required  for 
activation  of  rRNA  transcription  (Yuan  et  al.,  2007).  Thus,  we 
suggest  that  G9A-mediated  H3K9  methylation  serves  as  a  key 
link  between  amino  acid  sensing  and  transcriptional  control  of 
ribosome  biogenesis  and  cell  proliferation. 

Also  importantly,  we  show  that  G9A  requires  ATF4  for  tran¬ 
scriptional  activation  of  the  serine  biosynthetic  pathway  and  for 
stimulation  of  cell  proliferation.  ATF4  has  a  key  role  in  the  cellular 
response  to  amino  acid  limitation,  which  leads  to  increased 
ATF4  expression  primarily  at  the  translational  level.  ATF4  in 
turn  transcriptionally  activates  a  large  number  of  genes  including 


those  for  amino  acid  synthesis  and  trans¬ 
port  (Harding  et  al.,  2003;  Kilberg  et  al., 
2009;  Ye  et  al.,  201 2).  Although  the  mech¬ 
anistic  detail  remains  to  be  defined,  the 
observed  functional  connection  between 
G9A  and  ATF4  provides  further  evidence 
for  G9A  as  a  component  of  the  signaling 
pathway  that  responds  to  amino  acid 
abundance. 

It  is  well  documented  that  cancer  cells 
alter  cellular  metabolism  to  meet  the 
biosynthetic  challenge  of  growth  and  proliferation  (DeBerardinis 
et  al.,  2008).  Cancer  metabolism  is  characterized  by  aerobic 
glycolysis  with  a  high  rate  of  glucose  consumption  and  lactate 
production  (Warburg,  1956).  It  is  now  increasingly  recognized 
that  a  major  function  of  aerobic  glycolysis  is  to  divert  glycolic 
intermediates  for  the  biosynthesis  of  macromolecules  needed 
for  cell  growth  and  proliferation  (Cairns  et  al.,  2011;  Vander 
Heiden  et  al.,  2009).  Accumulated  evidence  suggests  that 
increased  activation  of  the  serine-glycine  biosynthetic  pathway 
is  an  important  part  of  cancer  metabolism  (Kalhan  and  Hanson, 
2012;  Locasale,  2013).  In  addition  to  generating  serine  and 
glycine  for  the  biosynthesis  of  proteins,  purines  (via  one-carbon 
units),  and  lipids  (via  phosphatidylserine),  this  pathway  pro¬ 
duces  equimolar  amounts  of  reduced  nicotinamide  adenine 
dinucleotide  (NADH),  a-KG,  and  5,10-MTHF  (Figure  7H).  These 
metabolites  have  critical  roles  in  the  control  of  cellular  meta¬ 
bolism  for  cell  proliferation  and  survival:  NADH  participates  in 
ATP  production  and  redox  regulation  (Corkey  and  Shirihai, 
201 2);  a-KG  supplies  carbon  to  the  tricarboxylic  acid  (TCA)  cycle 
for  the  generation  of  many  essential  biosynthetic  precursors 
(DeBerardinis  et  al.,  2008);  and  5,10-MTHF  is  a  coenzyme  for 
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Figure  7.  Activation  of  Serine-Glycine  Biosynthesis  Is  Required  for  G9A  to  Promote  Cell  Proliferation 

(A  and  B)  Immunoblotting  of  PHGDH  and  PSAT1  in  G9A-overexpressing  U20S  cells  with  PHGDH  (A)  or  PSAT1  (B)  knockdown.  PHGDH  and  PSAT1  levels 
were  quantified  against  a-tubulin. 

(C  and  D)  Growth  assay  of  G9A-overexpressing  U20S  cells  with  PHGDH  (C)  or  PSAT1  (D)  knockdown.  Error  bars  represent  SD  (n  =  3). 

(E  and  F)  Immunoblot  (E)  and  qRT-PCR  (F)  analysis  of  the  expression  of  serine-glycine  biosynthetic  enzyme  genes  in  G9A-overexpressing  U20S  cells  with 
ATF4  knockdown.  **p  <  0.001 . 

(G)  Growth  assay  of  G9A-overexpressing  U20S  cells  with  ATF4  knockdown.  Error  bars  represent  SD  (n  =  3). 

(H)  Model  for  G9A  regulation  of  serine  metabolism.  3-PG,  3-phosphoglycerate;  3-PHP,  3-phosphohydroxypyruvate;  3-PSer,  3-phosphoserine;  OAA, 
oxaloacetate;  Ptd-Ser,  phosphatidylserine. 

See  also  Figure  S5  and  Table  S3. 


the  only  cellular  pathway  of  de  novo  thymidylate  biosynthesis 
catalyzed  by  thymidylate  synthase  and  a  major  source  of  one- 
carbon  units  for  purine  synthesis  and  methyl  group  biogenesis, 
such  as  the  generation  of  SAM,  a  coenzyme  for  histone  and 
DNA  methyltransferases  (Kalhan  and  Hanson,  2012;  Teperino 
et  al.,  2010;  Tibbetts  and  Appling,  2010;  Touroutoglou  and 
Pazdur,  1996). 


Importantly,  we  present  evidence  suggesting  that  the  genera¬ 
tion  of  5,1 0-MTHF  is  crucial  for  the  function  of  the  serine-glycine 
biosynthetic  pathway  in  promoting  cancer  cell  survival  and 
proliferation.  We  show  that  supplemental  serine,  but  not  glycine, 
is  able  to  rescue  the  cell  death  phenotype  induced  by  G9A 
inhibition  or  silencing.  In  addition,  we  show  that  supplemental 
serine  synergizes  with  SHMT2  overexpression  to  enhance  not 
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only  cell  survival  but  also  cell  proliferation  under  the  condition  of 
G9A  inhibition.  Since  serine  and  glycine  are  interconvertible,  the 
most  obvious  explanation  for  these  observations  is  that  the 
conversion  of  serine  to  glycine  generates  5,10-MTHF,  a  metab¬ 
olite  essential  for  DNA  replication  and  one-carbon  metabolism. 
This  model  is  also  consistent  with  previous  reports  that  the 
mitochondrial  SHMT2  has  a  major  role  in  the  conversion  of  serine 
to  glycine  and  5,10-MTHF  (Herbig  et  al.,  2002;  Narkewicz  et  al., 
1 996;  Pfendner  and  Pizer,  1 980;  Stover  et  al.,  1 997;  Tibbetts  and 
Appling,  2010;  Yoshida  and  Kikuchi,  1970).  It  should  be  pointed 
out  that  the  mitochondrial  enzyme  glycine  cleavage  system 
catabolizes  glycine  into  carbon  dioxide,  ammonia,  and  5,10- 
MTHF  (Tibbetts  and  Appling,  2010).  Therefore,  supplemental 
glycine  could  be  a  source  of  5,10-MTHF  when  G9A  is  inhibited. 
However,  G9A  inhibition  markedly  reduces  the  steady-state 
levels  of  both  serine  and  glycine.  The  demand  for  serine  produc¬ 
tion  from  supplemental  glycine  would  decrease  the  5,10-MTHF 
pool.  Thus,  under  the  condition  of  G9A  inhibition,  5,10-MTHF 
production  from  the  cleavage  of  supplemental  glycine  might 
not  be  sufficient  to  sustain  DNA  synthesis  and  one-carbon 
metabolism. 

Recent  studies  have  identified  a  genetic  mechanism  for  acti¬ 
vating  the  serine-glycine  biosynthetic  pathway  in  breast  cancer 
and  melanoma  through  amplification  of  PHGDH  (Locasale 

et  al.,  2011;  Possemato  et  al.,  2011).  Our  findings,  coupled  to 
the  observation  that  G9A  is  overexpressed  in  various  cancers 

(Chen  et  al.,  2010;  Cho  et  al.,  2011;  Huang  et  al.,  2010;  Kondo 
et  al.,  2008),  suggest  an  epigenetic  mechanism  for  activating 
this  biosynthetic  pathway  in  cancer  cells.  Thus,  by  increasing 
the  amount  of  glycolic  carbon  diverted  into  macromolecule 
biosynthesis,  G9A  contributes  to  the  generation  of  a  meta¬ 
bolic  microenvironment  favorable  for  cancer  cell  survival  and 
proliferation. 

EXPERIMENTAL  PROCEDURES 
Cell  Culture  and  Reagents 

Cells  were  cultured  in  DMEM  (HyClone  SH30022,  Thermo  Scientific)  or  MEM 
(GIBCO  61 100-061,  Invitrogen)  supplemented  with  10%  dialyzed  FBS  (GIBCO 
26400-036).  BIX01294  (B931 1),  CQ  (C6628),  and  amino  acids  were  purchased 
from  Sigma-Aldrich.  Stock  solutions  of  10  mM  were  made  in  DMSO  (BIX)  or 
H20  (CQ).  In  all  BIX  experiments,  DMSO  was  used  as  control  (untreated). 
For  rescue  experiments,  amino  acids  (100  mM  stock  in  H20)  were  added  to 
MEM  at  the  final  concentration  of  0.4  mM  (same  as  in  DMEM).  For  serine  depri¬ 
vation  assay,  cells  were  cultured  in  MEM  supplemented  with  nonessential 
amino  acids  (MEM-NEAA)  and  then  transferred  to  serine-deficient  MEM- 
NEAA.  Phase  contrast  images  were  captured  using  an  Axio  Observer 
microscope  and  AxioVision  software  (Carl  Zeiss  Microimaging).  Trypan  blue 
exclusion  assay  was  used  to  assess  cell  growth  and  survival. 

Immunoblotting 

Cells  were  suspended  in  SDS  sample  buffer  and  protein  concentrations 
determined  using  a  protein  assay  kit  (Bio-Rad).  Proteins  (20  ^g)  were  sepa¬ 
rated  on  SDS-polyacrylamide  gels,  transferred  to  nitrocellulose  membranes, 
and  probed  with  primary  antibodies  detailed  in  Supplemental  Experimental 
Procedures.  Horseradish  peroxidase-conjugated  goat  anti-mouse  and  anti¬ 
rabbit  IgG  (Santa  Cruz  Biotechnology)  were  used  as  secondary  antibodies. 
Proteins  were  visualized  using  a  SuperSignal  West  Pico  chemiluminescence 
kit  (Pierce)  and  quantified  with  ImageJ  (version  1.47d).  Films  were  exposed 
for  various  times  for  protein  quantification  within  linear  range  of  detection. 
For  visualization  and  quantification  with  the  Odyssey  system,  goat  anti-mouse 
IRDye  800  or  680  and  anti-rabbit  IRDye  800  or  680  (LI-COR  Biosciences)  were 
used  as  secondary  antibodies. 


Immunofluorescence 

Cells  were  fixed  with  4%  paraformaldehyde  and  costained  with  rabbit  anti- 
LC3B  (3868,  Cell  Signaling,  1:200)  and  mouse  anti-a-tubulin  (B-5-1-2, 
Sigma-Aldrich,  1:2000).  All  secondary  fluorescence  antibodies  (Alexa  Fluor 
488  and  594)  were  from  Molecular  Probes  and  used  at  1 :1 ,000  dilution.  Nuclei 
were  stained  with  DAPI  (D3571,  Invitrogen)  or  Hoechst  33342  (H3570, 
Invitrogen).  Fluorescent  images  were  captured  with  an  Axio  Observer 
microscope  and  AxioVision  software.  For  quantification  of  LC3B-positive 
puncta  per  cell,  approximately  100  cells  were  counted  from  a  random 
selection  of  four  to  six  630  x  fields  independently  by  two  investigators. 

Microarray 

Total  RNA  was  isolated  using  Trizol  (Invitrogen)  from  three  independent  sam¬ 
ples  of  BE(2)-C  cells  with  or  without  5  [xM  BIX  for  24  hr.  Affymetrix  microarray 
was  performed  using  the  Human  Gene  1.0  ST  microarray  chip.  Data  were 
normalized,  significance  determined  by  ANOVA,  and  fold  change  calculated 
with  the  Partek  Genomics  Suite.  GO  analysis  was  performed  with  DAVID 
(Huang  et  al.,  2009)  for  all  differentially  expressed  genes  (>  ±1 .5-fold,  p  <  0.01). 

RNA-Seq 

Total  RNA  was  isolated  using  Trizol  from  three  independent  samples  of 
U20S  cells  either  untreated  or  treated  with  5  jxM  BIX  for  6  or  24  hr.  RNA-seq 
libraries  were  generated  with  an  lllumina  TruSeq  RNA  sample  preparation  kit 
(RS-1 22-2001)  and  sequenced  using  lllumina  high-seq  2000  with  a  read  length 
of  50  bp  with  pair  ends.  RNA-seq  reads  were  mapped  to  the  human  genome 
(hg19)  using  TopHat  (Trapnell  etal.,  2009).  Only  those  reads  mapped  to  unique 
genomic  locations  and  with  <5%  mismatches  were  analyzed  further.  We  used 
Cufflinks  to  measure  gene  transcripts,  and  Cuffdiff  to  identify  differentially 
expressed  genes  with  a  RefSeq  GTF  file  downloaded  from  lllumina  iGenomes 
(Roberts  et  al.,  201 1).  Deferentially  expressed  genes  were  fed  into  DAVID  for 
functional  annotation.  Bioconductor  (Gentleman  et  al.,  2004)  packages  were 
used  to  generate  heatmap. 

qRT-PCR 

Total  RNA  was  isolated  from  three  independent  samples  using  Trizol.  Reverse 
transcription  was  performed  using  Superscript  II  Reverse  Transcriptase 
(Invitrogen).  Quantitative  real-time  PCR  was  performed  using  a  RT2  SYBR 
green/Fluorescein  PCR  master  mix  (SABiosciences)  on  an  iQ5  real-time 
PCR  system  (Bio-Rad)  with  primers  against  ATF4,  CCNA2,  CCNB1 ,  CCNB2, 
CDC25C,  PHGDH,  PSAT1 ,  PSPH,  SHMT1 ,  SHMT2,  and  B2M  (Table  S3).  All 
primer  pairs  were  verified  by  melting  curve  analysis  following  qRT-PCR,  with 
each  primer  pair  showing  a  single  desired  amplification  peak.  All  samples 
were  normalized  to  (32  microglobulin  mRNA  levels. 

ChIP-qPCR  and  ChIP-Seq 

ChIP  was  performed  as  described  (Lee  et  al.,  2006),  using  2  x  107  parental, 
GFP-overexpressing,  or  G9A-overexpressing  U20S  cells  with  or  without 
5  [xM  BIX  for  the  indicated  times.  Crosslinked  chromatin  was  sheared  through 
sonication  and  immunoprecipitated  using  ChIP  grade  mouse  anti-H3K9me1 
(17-680),  rabbit  anti-H3K9me2  (17-648),  or  control  mouse  IgG  or  rabbit  IgG 
(all  from  Millipore).  For  qPCR,  two  independent  ChIP  samples  were  analyzed, 
and  each  sample  was  assayed  in  triplicate  using  primers  that  cover  the  pro¬ 
moter  regions  of  PHGDH  and  PS ATI ,  and  the  coding  region  of  GAPDH 
(detailed  in  Supplemental  Experimental  Procedures).  For  ChIP-seq,  libraries 
were  generated  from  ChIP  samples  using  an  lllumina  TruSeq  ChIP  Sample 
Prep  kit  (IP-202-1012)  and  sequenced  using  lllumina  high-seq  2000  with  a 
read  length  of  50  bp  with  pair  ends.  Raw  lllumina  sequencing  reads  in  the 
FASTQ  format  were  cleaned  using  in-house  scripts  by  trimming  sequencing 
adaptors  and  low-quality  bases  in  both  ends  (Q  <  67  in  lllumina  1 .5).  Cleaned 
sequences  were  then  mapped  to  the  human  genome  (hg19)  using  Novoalign 
v2.07.  The  reads  that  mapped  uniquely  to  a  single  genomic  locus  were  used 
for  peak  finding  with  MACS  vl.4  (Zhang  et  al.,  2008),  and  only  those  peaks 
with  false  discovery  rate  (FDR)  <1  %  were  compared  with  RefSeq  genes  in 
the  UCSC  genome  browser. 

Metabolite  Analysis 

GC-MS  metabolomic  analysis  was  performed  in  the  University  of  Utah 
Metabolomics  Core.  Parental,  GFP-  or  G9A-overexpressing  U20S  cells 
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were  cultured  in  MEM  supplemented  with  10%  dialyzed  FBS  and  were  either 
untreated  or  treated  with  5  ^M  BIX  for  4  hr.  Cells  were  washed  with  PBS,  har¬ 
vested  by  trypsin  digestion,  transferred  to  a  microfuge  tube,  and  frozen  on  dry 
ice.  Metabolites  were  extracted  by  80%  methanol  at  -20°C  and  dried  by 
vacuum  centrifugation.  Six  biological  replicate  samples  (5  x  106  cells/sample) 
were  analyzed  for  each  condition.  GC-MS  analysis  was  performed  with  a 
Waters  GCT  Premier  mass  spectrometer  fitted  with  an  Agilent  6890  gas 
chromatograph  and  a  Gerstel  MPS2  autosampler.  Data  were  collected  using 
MassLynx  4.1  software  (Waters).  Metabolites  were  identified  and  their  peak 
area  was  recorded  using  QuanLynx.  Data  were  normalized  for  extraction  effi¬ 
ciency  and  analytical  variation  by  mean  centering  the  area  of  D4-succinate. 

Flux  Analysis 

U20S  cells  with  or  without  overexpression  of  GFP  or  G9A  were  cultured  in 
either  glucose-deficient  DMEM  (11966-025,  Invitogen)  or  MEM  supplemented 
with  10%  dialyzed  FBS  and  25  mM  [U-13C]  glucose  (CDLM-3813-1 ,  Cam¬ 
bridge  Isotope  Laboratories)  with  or  without  5  fiM  BIX  for  16  hr  before  metab¬ 
olite  extraction  with  80%  methanol.  Conversion  of  [U-13C]  glucose  to  other 
metabolites  was  analyzed  by  targeted  LC-MS/MS  using  selected  reaction 
monitoring  (SRM)  with  a  5500  QTRAP  mas  spectrometer  (Yuan  et  al.,  2012). 
Three  biological  replicate  samples  were  analyzed  for  each  condition. 

Analysis  of  Tumor  Data  Sets 

The  Versteeg  NB88  data  set  contains  88  primary  neuroblastoma  tumors  of  all 
stages  (Valentijn  et  al.,  2012).  Kaplan-Meier  analysis  was  conducted  online 
(http://r2.amc.nl),  and  the  resulting  survival  curve  and  p  value  (log-rank  test) 
were  downloaded. 

Statistics 

Data  are  presented  as  mean  ±  SD.  Statistics  were  determined  with  unpaired, 
two-tailed  Student’s  t  test  using  GraphPad  Prism  6.01  for  Mac. 

ACCESSION  NUMBERS 

The  NCBI  Gene  Expression  Omnibus  (GEO)  accession  number  for  the  micro¬ 
array  data  reported  in  this  paper  is  GSE51512. 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  five  figures,  three  tables,  and  Supple¬ 
mental  Experimental  Procedures  and  can  be  found  with  this  article  at  http:// 
dx.doi.org/1 0. 1 01 6/j.cmet.201 3.1 1 .004. 
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